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Preface 


The  purpose  of  Che  thesis  was  Co  develop  a  model  for  Che  forced 
transverse  response  of  a  misCuned  disk  chat  could  predict  rotating  node 
lines  near  a  resonance.  A  stationary  disk  was  used  because  the  rotating 
disk  equation  of  motion  has  not  been  solved.  The  disk  mistuning  was 
Introduced  by  placing  torsional  springs  at  the  disk  outer  edge  but  only 
at  the  nodes  of  given  modes.  This  selective  location  of  the  springs 
keeps  the  analytic  representation  of  the  disk's  response  very  simple. 
Then,  the  boundary  conditions  were  specified  as  having  a  spring 
resistance  only  at  these  points,  and  another  set  of  boundary  conditions 
were  specified  as  free  only  at  the  points  of  maximum  displacement  on  the 
edge .  This  rather  contrived  form  of  imperfection  was  needed  to  keep  the 
eigenfunctions  simple  enough  for  a  closed  form  solution  to  the  response. 
Thus,  our  mistuning  model  is  not  at  all  a  true  analytical  representation 
of  an  imperfect  disk.  It  is  really  just  a  mathematical  trick  used  to 
introduce  an  asymmetiry  into  the  boundary  conditions,  and  it  is  not 
expected  to  be  useful  in  the  modeling  of  any  particular  real 
Imperfection.  What  it  does,  however,  is  allow  a  qualitative 


understanding  of  how  a  real  mistuned  disk  might  respond  to  certain 
loads.  In  this  respect,  it  is  better  than  a  true  model  of  an  imperfect 
disk.  Because  the  mistuned  mode  shapes  observed  in  experiments  look 
just  like  those  of  the  eigenfunctions  calculated  in  this  thesis,  I 
believe  that  most  forms  of  real  mistuning  will  have  responses  similar  to 
the  closed  form  results  given  herein  (see  1:12) 


Our  mistuned  disk  response  model  cannot  evaluate  the  response  of  a  _ 

n  For 

design.  The  actual  mistuned  resonant  frequencies  of  the  disk  as  well  as  - 

the  modal  fractions  of  critical  damping  must  be  provided  for  the  model 

ed 

to  be  useful.  Thus,  the  actual  disk  must  exist,  and  experiments  must  be  tiQn_ 
performed  to  obtain  these  modal  data  before  the  model  is  applied. 
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The  characteristics  of  the  mistimed  stationary  disk  response  also 
appear  to  be  attributable  to  a  mistimed  rotating  disk.  This  is  because 
the  mode  shapes  of  the  rotating  disk  are  assumed  to  have  a  sinusoidal  9 
dependence,  which  is  the  same  form  as  that  of  the  stationary  disk,  and 
the  stationary  disk  response  model  has  the  same  characteristics  for  all 
such  mode  shapes. 

I  would  like  to  thank  my  wife  for  her  patience  and  support  during 
this  study.  I  would  also  like  to  thank  Lt  Col  Bagley,  one  of  the  best 
thesis  advisors  at  AFIT,  for  his  expertise,  insight,  and  comprehensive 
support.  I  am  also  indebted  to  Lt  Col  Baker  of  the  Math  Department  for 
his  help  in  proving  the  self-adjointedness  of  the  structural  models  used 
in  this  thesis. 
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Abstract 


The  objective  of  this  thesis  was  to  develop  a  simple  model  for  the 
forced  response  of  a  mistuned,  center- clamped  disk  to  help  predict 
previously  unexplained  experimental  observations.  Experiments  have 
shown  that  the  response  near  a  resonance  can  take  the  form  of  rotating 
node  lines  (1:12).  A  simple  model  of  a  mistuned  stationary  disk  was 
developed  and  the  transverse  responses  to  both  rotating  point  loads  and 
stationary  point  loads  oscillating  in  time  were  determined.  The  disk 
model  uses  torsional  springs  applied  at  points  on  the  outer  edge  of  a 
stationary  annular  disk  to  mistune  a  given  transverse  vibration  mode 
such  that  the  two  resulting  mistuned  mode  shapes  are  each  described  by 
only  four  terms.  The  response  of  all  modes  other  than  the  two  mistuned 
modes  being  examined  was  neglected,  as  only  the  behavior  near  a 
resonance  was  considered. 

The  responses  to  both  load  types  did  indeed  have  rotating  node 
lines  near  a  resonance.  The  node  lines  were  predicted  to  be  slightly 
distorted  and  rotating  at  a  varying  but  average  angular  speed  identical 
to  the  circular  frequency  of  excitation.  In  addition,  it  was  predicted 
that,  if  the  modal  response  amplitudes  are  equal  and  the  modal  responses 
are  ir/2  radians  out  of  phase,  a  node  point  rotating  at  the  constant 
excitation  speed  would  occur.  The  results  also  indicate  that  the  node 
lines  can  rotate  in  either  direction  due  to  a  stationary  load  but  only 
backward  due  to  a  fot^rard  rotating  load. 

In  a  qualitative  comparison,  the  resp'^nse  of  a  rotating  disk  to  a 
stationary  point  load  was  found  to  be  similar  to  the  stationary  disk 
response  to  a  rotating  load. 
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A  NEW  MODEL  FOR  THE  FORCED  RESPONSE  OF  A  MISTUNED  DISK 


I.  Introduction 

A  perfect  disk  vibrating  transversely  can  theoretically  have 
diametrical  node  lines  moving  through  it  or  remaining  stationary  at  an 
arbitrary  location,  but  a  real  disk  can  only  have  stationary  node  lines 
or  node  lines  moving  at  Integer  fractions  of  the  excitation  speed.  The 
difference  is  due  to  imperfections  in  the  real  disk  that  split  any  given 
mode  into  two  "mlstuned"  modes  that  each  have  their  own  stationary  node 
lines.  The  reason  the  node  lines  of  the  real  disk  are  able  to  move  is 
that  the  disk  response  includes  the  response  of  both  mlstuned  modes, 
which  may  be  out  of  phase . 

in  Free  Vibration 

The  stationary  disk  equation  of  motion  and  its  solution  are  well 
known;  however  such  is  not  the  case  with  the  rotating  disk.  Southwell 
published  the  natural  frequencies  of  the  certer-clamped  annular  plate  in 
1921  (2:139).  In  1970,  Mote  published  the  eigenfunction  coefficients 
for  the  center -clamped  plate  (3:332-333).  The  rotating  disk  equation  of 
motion  was  studied  by  Lamb  and  Southwell  in  1921  (4:276-280).  They  were 
able  to  estimate  the  natural  frequencies  of  the  rotating  disk  but  could 
not  solve  the  equation  of  motion  in  closed  form.  In  1972,  Barasch  and 
Chen  verified  the  Lamb  and  Southwell  estimates  using  numerical 
techniques;  however,  no  one  has  yet  determined  the  closed  form  solution 
(5:1143) . 

Research  in  Forced  Response 

In  1957,  Tobias  and  Arnold  postulated  that  a  node  line  could  move 
through  a  perfect  stationary  disk  at  the  speed  d9/dt  -  w/n  ,  where  w  is 
the  excitation  frequency  and  n  is  the  number  of  nodal  diameters  in  the 
given  mode;  however,  they  made  no  such  prediction  for  the  mistuned  disk 
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(6:672).  Then,  in  1981,  Stange  and  MacBain  reported  the  observation  of 
such  a  response  In  a  mlstuned  disk  (1:5).  This  observation  led  to 
further  study  and,  eventually,  to  this  thesis. 

One  experiment  of  Stange  and  MacBain  involved  exciting  a 
stationary,  center-clamped,  bladed  disk  with  a  stationary  electromagnet¬ 
ic  exciter.  They  excited  the  modes  with  two  nodal  diameters  and  no 
nodal  circles.  The  exciter  was  located  midway  between  adjacent  node 
lines  of  the  two  mlstuned  modes  so  as  to  excite  both  modes  equally.  In 
describing  the  response  due  to  excitation  between  the  two  mistimed 
resonant  frequencies,  they  wrote  the  following: 

By  exciting  the  disk  at  a  frequency  between  these  two 
frequencies  though,  a  point  could  be  found  at  which  the 
dominant  strain  gages  for  each  mode  became  equal  in  amplitude. 

Double  pulsed  holograms  taken  at  this  point  consistently  showed 
a  2N  mode,  but  in  different  orientations  with  respect  to  the 
disk,  Indicating  that  the  mode  was  rotating.  . . .  Switching  the 
laser  into  its  strobed  mode,  it  was  found  that  by  strobing  the 
disk  at  exactly  the  excitation  frequency,  the  rotating  2N  mode 
could  be  made  to  appear  stationary.  Changing  the  phase  of  the 
strobe  resulted  in  a  change  in  the  orientation  of  the  mode, 
thus  showing  that  the  mode  was  rotating  with  respect  to  the 
disk  at  the  excitation  frequency.  (1:5) 

The  mlstunlng  in  this  experiment  was  due  only  to  the  manufacturing 
imperfections  (it  was  not  induced).  These  imperfections  were  small  as 
measured  by  the  amount  of  mistuning,  given  by 

(,u>2  -  Wj)/w^  -  (266.5  Hz  -  263.6  Hz)/263.6  Hz  -  0.011  .  The  researchers 
did  not  report  any  knowledge  of  the  exact  form  of  imperfection  that 
caused  the  mistuning,  and  it  would  be  surprising  if  anyone  could. 

Because  of  the  mathematical  difficulties  involved  with  imperfect 
disks,  most  of  the  mlstuned  disk  theoretical  research  has  been 
numerical.  Imperfections  tend  to  require  infinite  series  representa¬ 
tions  of  the  eigenfunctions.  Numerical  methods  do  not,  however,  yield 
much  insight.  Thus,  an  analytical  approach  is  preferred  in  searching 
for  the  cause  of  the  rotating  node  lines.  An  analytical  approach, 
however,  will  require  some  care  in  the  modeling  of  the  mistuning  and  may 
require  a  departure  from  the  true  modeling  of  a  particular  t3rpe  of 
imperfection. 
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Available  Migtwlng  Model 

Marita  and  Lelssa  have  shown  that  mlstuning  can  be  modeled  by 
introducing  non-uniform  boundary  conditions  (7:109).  In  the  referenced 
work,  the  authors  apply  rotatory  springs  to  the  simply- supported  outer 
boundary  of  a  disk.  This  approach  leads  to  an  infinite  series  of  Bessel 
functions  representation  of  the  mode  shapes.  Furthermore,  it  is 
apparent  from  the  referenced  work  that  such  coupling  is  almost 
inevitable  when  the  mlstuning  is  introduced  through  the  arbitrary 
addition  of  external  mass  or  stiffness.  This  is  indicated  by  the 
movement  of  the  mode  node  line  off  the  disk  center  for  one  of  the  split 
modes.  Figure  1  is  a  sketch  of  one  result  given  by  Marita  and  Leissa 
(7:108). 


Figure  1.  Asymmetric  Node  Line  of  a  Mistuned  Disk 
Mode:  m  -  0,  n  -  1  (Taken  from  7:108) 


Such  a  node  line  cannot  exist  in  a  mode  shape  of  the  form 
f(r)  sin  n9  because  the  angular  dependence  of  the  displacement  is 
clearly  not  sinusoidal.  In  fact,  the  only  node  forms  possible  In  the 
above  solution  are  nodal  diameters  and  nodal  circles.  Thus,  a  more 
refined  mlstunlng  model  Is  needed. 

A  desirable  model  for  the  purpose  of  this  thesis  would  be  one  that 
Introduces  the  Imperfection  without  coupling  the  Bessel  functions.  Then 
the  mode  shapes  could  be  described  by  a  few  terms,  and  the  response  of 
the  two  mlstuned  modes  may  possibly  be  added  without  concern  about  all 
other  orders  of  the  Bessel  functions.  In  light  of  the  previous 
research,  the  obvious  place  to  Introduce  Che  Imperfection  Is  at  Che  disk 
boundary . 

The  type  of  Imperfection  must  also  be  carefully  considered.  If  a 
mass  Imperfection  Is  added  at  the  boundary,  Che  model  may  be  more 
accurate;  however,  the  boundary  condition  will  Involve  the  eigenvalue, 
and  the  resulting  eigenvalue  problem  will  not  be  self-adjoint.  As  a 
direct  result,  Che  mode  shapes  become  Infinite  series  of  eigenfunctions. 
A  stiffness  Imperfection  Is  better  In  this  sense,  since  springs  at  the 
boundary  do  not  usually  result  In  an  eigenvalue  dependent  boundary 
condition. 
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II.  Stationary  Disk  Equation  2f  fiDSi  Solutions 


The  equation  of  notion  and  corresponding  boundary  conditions  of 
the  stationary  annular  plate  are  developed  here  using  Hamilton's 
Principle.  The  development  neglects  the  effects  of  transverse  shear  and 
rotatory  Inertia,  and  Is  therefore  limited  In  application  to  thin  plates 
and  the  lower  few  resonant  modes  of  vibration. 

Motion 

The  plate  strain  and  kinetic  energy  for  the  above  mentioned 
restrictions  are  given  by  Uarburton  (8:242).  They  are 
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E  -  Young's  Modulus 
h  -  the  disk  thickness 
u  •  Poisson's  Ratio 
p  -  the  disk  volumetric  density 

u;  -  the  disk  transverse  deflection 

a  •  the  disk  Inner  radius 
h  -  the  disk  outer  radius 
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The  work  done  on  the  plate  by  a  transverse  load  q  is 
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Hamilton's  Principle  states  that 
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where  the  S  refers  to  the  variation. 


(3) 


Substituting  equations  1,  2,  and  3  into  equation  4  and  performing  the 
variation  yields  the  equation  of  motion  for  the  forced  response  of  a 
stationary  disk, 
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and  the  stationary  disk  boundary  conditions 
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r^dG^*  dG 


0 


(8) 


L 

D 

d^W  d^w  \  dw  (2-v)  d^w 

{v-3)d^w' 

6w 

dr^  dr^  r  dr  r  drdG^ 

r^  dG^  _ 

(9) 


L 


6 


2£>(1  -V) 


[r^3e~r3rc 


These  boundary  conditions  require  sone  interpretation.  The 
conditions  expressed  by  equations  7  and  8  apply  to  a  boundary  where  9  is 
constant.  In  the  case  of  an  annular  disk,  these  equations  are  merely  an 
expression  of  the  continuity  of  the  disk  in  the  9  direction  and  do  not 
require  application.  The  condition  expressed  by  equation  10  applies 
only  to  comers  and  is  usually  satisfied  by  solutions  consistent  with 
the  physical  boundary  conditions  given  by  equations  6  and  9.  The 
remaining  eqviatlons  (6  and  9)  are  the  moment  and  shear  conditions  on  the 
inner  and  outer  circumferences.  Equation  6  states  that  either  the  slope 
of  the  transverse  deflection  or  the  moment  must  be  specified  at  each 
circumferential  boundary.  Equation  9  states  that  either  the  shear  or 
the  transverse  deflection  must  be  specified  at  each  circumferential 
boundary.  After  selecting  a  set  of  boundary  conditions,  these  equations 
must  be  applied  to  the  general  solution  of  the  equation  of  motion  to 
determine  the  specific  solution  to  a  given  disk  and  its  supports. 


SolMtl9n  o£  Equation  Motion 

To  determine  the  forced  response  of  the  continuous  structure,  the 
solutions  to  the  homogeneous  equation  must  be  found  first.  This  is  done 
using  the  method  of  separation  of  variables.  Then,  the  solution  to  the 
non-homogeneous  equation  is  assumed  to  be  the  sum  of  all  of  the 
solutions  to  the  homogeneous  equation.  Meirovitch  calls  this  the 
expansion  theorem  (9:143). 

The  homogeneous  form  of  the  stationary  disk  equation  of  motion  is 
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Assuming  the  time  and  spatial  variables  can  be  separated,  let 

wir.e.t)^W(r.e)aCt)  (12) 

Substituting  this  equation  into  eq  11  and  separating  variables  yields 

V^l2-^co^i2J-0  (13) 

Letting  phu/^/D  -  yields 


(v"-^")«Z(r.e)-0 


or 


(14) 


The  solution  to  this  equation  is 


*  D  „^r]]s\nne  *[E  F  G 

„^K^[p^^r)]cosn9  d^) 

where 

-  the  mode  eigenvalue 

J„  -  the  Bessel  function  of  the  first  kind  of  order  n 

-  the  Bessel  function  of  the  second  kind  of  order  n 

I„  -  the  modified  Bessel  function  of  the  first  kind  of  order  n 

-  the  modified  Bessel  function  of  the  second  kind  of  order  n 

-  the  eigenfunction  coefficients  to  be  determined  by  applying 
the  boundary  and  initial  conditions 

and  the  mn  mode  has  m  nodal  circles  and  n  nodal  diameters. 
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III.  Mis tuning  Model 


To  obtain  a  simple  expression  for  the  stationary  disk  response,  a 
simple  mistvmlng  model  must  be  used.  Unfortxinately ,  as  mentioned 
previously,  such  a  model  Is  unavailable  In  the  literature.  A  simple 
model  can  be  found  If  the  boundary  conditions  used  to  Introduce  the 
Imperfection  are  chosen  carefully.  Care  should  also  be  taken  to  ensure 
the  orthogonality  of  the  mistvined  modes  so  that  the  differential 
equations  of  motion  will  not  be  coupled.  The  approach  will  be  to  apply 
springs  to  the  boundary  such  that  eight  boundary  condition  equations 
result  and  yield  two  simply  expressed  mlstuned  modes  for  each  mode  of  a 
perfectly  symmetric  disk. 

Simple  Model 

Consider  the  annular  disk  shown  In  figure  2.  The  disk  Is  clamped 
on  the  Inner  boundary  (r  -  a)  with  torsional  springs  applied  at  certain 
points  on  the  outer  boundary  (r  -  b) .  Let  the  outer  boundary  be 
specified  as  free  at  certain  other  points  (b.Sf),  where  Jp  -  it/2n  +  iit/n 
(i  -  0,l,2,...,2n  •  1),  but  unspecified  at  all  remaining  points.  Let 
the  springs  be  of  equal  stiffness  K  and  be  applied  at  all  points  (b,^^), 
where  tfg  -  iw/n  (i  -  0,1,2, ...,2n  -  1).  The  number  of  springs  required 
is  obviously  dependent  on  the  number  of  nodal  diameters  n  in  the  mode 
being  mlstuned.  Figure  2  Illustrates  the  case  n  -  2.  Also  consider  in 
this  figure  that  the  torsional  spring  mounts  are  massless  and  are  free 
to  move  transversely  with  the  disk  but  not  free  to  rotate  with  the  disk. 

The  boundary  conditions  applied  In  this  model  will  mlstune  the 
disk  because  they  create  two  different  boundary  condition  equations  for 
the  outer  disk  edge.  One  mode,  having  its  node  lines  at  the  angles  6^ 
will  have  to  rotate  against  the  torsional  springs,  since  they  are 
located  at  its  lines  of  peak  displacement.  Another  mode  will  exist  but 
will  not  be  affected  by  the  springs  because  they  are  attached  at  its 
node  lines. 


Figure  2.  Clanped  Annular  Disk  With  Torsional  Springs 
Applicable  to  Modes  With  n  -  2 


The  boundary  conditions  are  described  by  the  following  eight  equations: 


jr(a  ,0j)  -  0 


(16) 


dOj/  ^ 

—  (a.ej-o 


(17) 


A/(b.0,)-A:^(b.0,) 


(18) 
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(19) 

a}(a.0f)-O 

(20) 

diD (  .  \  - 

—  (a,0,)-O 

(21) 

A/(b.0,)-O 

(22) 

V'(b.fl,)-0 

(23) 

where 

M  -  the  moment  applied  to  the  boundary  (eq  6) 

V  -  the  transverse  shear  applied  to  the  boundary  (eq  9) 

K  -  the  torsional  spring  constant 

Notice  that  the  above  equations  give  the  same  result  for  all  values  of 
the  index  i  because  each  of  these  equations  is  of  the  form 
fi(r)  sin  n$  +  /^(r)  cos  n$  -  0  and  cos  ntfg  -  ±1,  sin  n$f  -  ±1, 
cos  nSf  •  0,  and  sin  ntfg  -  0  for  all  values  of  I.  Observe  also  that 
the  boundary  condition  points  of  application  are  chosen  so  as  to  create 
two  independent  sets  of  four  equations.  This  can  be  seen  by  considering 
eqxiation  15  and  observing  that  cos  ntfp  -  0  and  sin  n^g  -  0  .  Thus, 
equations  20-23  can  be  used  to  solve  for  the  first  four  coefficients  in 
equation  15,  and  equations  16-19  can  be  used  to  solve  for  the  last  four 
coefficients.  This  will  result  in  two  orthogonal  modes  of  a  single 
order  n,  each  described  by  only  four  Bessel  functions.  Thus,  these 
boundary  conditions  have  been  specifically  chosen  to  produce  separate 
mode  shapes  that  need  only  a  few  Bessel  functions  to  describe  them. 

Stating  the  boundary  conditions  in  this  fashion  leaves  no  means  of 
mathematically  expressing  the  condition  on  the  outer  boundary  except  at 
the  points  (b.fip)  and  (b,fg);  however,  it  provides  a  set  of  eight 
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equations  that  can  be  used  to  solve  the  eigenvalue  problem. 

Selecting  the  boundairy  conditions  from  eqxiatlons  6  and  9  and 
sxibstltutlng  them  Into  equations  16-23  and  simplifying  results  in  the 
following  eight  equations: 


COS/I0J.-O  (24) 

■*‘^AtA|^^  ~  ^  nn^)  ~  ^  n*  I  (  fi  nn^  ]  | 

cosa0j  =  O  (25) 


mn 


mn^) 


*f™{'>y.(?-.t>)-P„t>K..,(^„6) 

P  mu  ^ 


~  ^  nni^^  n{fi  nn^)'*'  P  »*  l{0  nn^] 

~  ^  nni^^  ni^  mn^]^  ^  n*\{fi  n»^) 


( 1  - v)n^ 


[(A-  0/f«(^„.b)-^„.by,.,(^^,b)]}]cos;ie^-0  (27) 


svnn0,«O  (28) 


\_H  nix'*  J  |_p„„a 

L^'mn^  J  LP^„U 


sinn0-  =  O  (29) 


sinn0f“O  (^0) 

[/l..l'l-'.(^».6)-#~<>-'..,(/S™6) 

P  mu  O 

*B„.{ny,(fi^b)-0^br,.,{fi^b] 

P  mn^ 

-C^{nl,(fi^b)*0^bl,.,(fi^b) 

-^^^[('•-l)/.(P..b)-?„.b;,.,(p^b)]} 

H  «n  ^ 

-D^.{nK.[p„,b)-$^bK..,(e^b) 

P  mn  ^ 


Equations  15  and  2^-31  are  used  to  construct  the  solution  for  the 
free  vibration  of  the  mistuned  plate  shown  in  figure  2.  They  are 
precisely  the  equations  for  two  separate  structures,  namely  the 
clamped- sprung  annular  plate  (represented  by  eqs  24-27)  and  the 
clamped-free  annular  plate  (represented  by  eqs  28-31);  and,  they  create 
two  distinct  modes,  because  the  boundary  conditions  of  one  mode  are 
applied  at  the  nodes  of  the  other.  That  is,  the  two  sets  of  points 
where  the  boundary  conditions  are  applied  are  nodes  of  the  other 
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mistimed  disk  mode  shape.  The  specification  of  the  outer  boundary 
condition  only  at  these  points  is  what  allows  the  two  different  mode 
shapes  of  the  disk  to  be  independent  of  each  other.  Equations  28-31  are 
similar  to  those  developed  by  Southwell  (2:138-139). 

As  mentioned  previously,  since  the  eight  characteristic  equations 
can  be  separated  into  two  independent  sets  of  four  characteristic 
equations,  the  eigenvalue  and  eigenfunction  coefficients  for  each  set  of 
equations  can  be  solved  for  independent  of  the  other  eigenvalue  and 
eigenfunction  coefficients.  The  two  resulting  eigenvalues  will 
converge  to  a  single  eigenvalue  in  the  limit  of  zero  spring  stiffness, 
as  is  obvious  when  one  compares  equation  26  to  equation  30  (note  that 
the  other  three  equations  in  each  set  of  four  are  identical).  Thus,  the 
degree  of  mistunlng  in  the  model  is  directly  related  to  the  spring 
stiffness,  and  a  zero  spring  stiffness  corresponds  to  a  perfect  disk. 

Henceforth,  the  mode  produced  by  the  clamped- free  boundary 
conditions  will  be  designated  mode  1,  or  the  sine  mode.  Alternately, 
the  mode  produced  by  the  clamped- sprung  plate  will  be  designated  mode  2, 
or  the  cosine  mode.  The  two  mode  shapes  are  represented  as 

and 

Split  Modes  in  £  Tuned  DisJc 

In  the  modeling  of  the  mistiming  phenomenon,  it  is  desirable  to 
show  that  the  mistunlng  can  be  reduced  as  the  amount  of  imperfection  is 
reduced.  This  can  be  shown  by  reducing  the  degree  of  imperfection  in 
the  structure  and  observing  the  corresponding  degree  of  difference 
between  the  natural  frequencies  and  mode  shapes  of  the  two  mistuned 
modes.  In  the  model  developed  in  the  previous  section,  the  degree  of 
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imperfection  is  a  function  of  the  spring  constant  K.  An  example  problem 
will  illustrate  the  convergence  of  the  two  mistuned  modes  in  the  limit 
of  small  K. 

Equations  IS  and  24- 31  can  be  solved  to  yield  the  two  mistuned 
mode  shapes  and  frequencies  for  each  set  of  values  m,  n.  To  solve  for 
the  eigenvalues  and  eigenfunction  coefficients,  the  matrix  of  the 
coefficients  of  the  boundary  condition  equations  must  be  formed.  This 
matrix  consists  of  terms  such  as  a^i  -  (/3„a)  and  a^j  -  (/Jan®)  • 

The  matrix  is  a  4X4  for  each  of  the  two  modes .  Note  that  only  one  row 
will  be  different  in  each  of  the  two  4X4  matrices  since  only  one 
boundary  condition  is  different.  The  eigenvalues  for  given  values  of  a, 
b,  and  n  are  determined  first.  This  is  done  by  estimating  an  eigenvalue 
and  calculating  the  determinant.  If  the  determinant  is  not  zero,  then 
the  eigenvalue  is  incremented  and  the  determinant  re-calculated  until 
found  to  be  approximately  zero.  The  lowest  eigenvalue  for  which  the 
determinant  is  zero  corresponds  to  the  m  -  0  mode,  and  each  successive 
zero  of  the  determinant  corresponds  to  a  successively  higher  value  of  m. 

Once  the  eigenvalues  are  determined,  the  coefficients  of  the 
eigenfunctions  can  be  calculated.  Since  one  coefficient  will  be 
arbitrary,  the  first  coefficient  may  be  assigned  the  value  of  1.  This 
will  leave  only  three  unknowns  for  each  set  of  four  equations.  The 
remaining  coefficients  may  be  determined  by  substituting  the  appropriate 
eigenvalue  into  any  three  of  the  four  equations  and  using  any  equation 
solving  technique.  This  approach  leads  to  sine  mode  coefficients 
^nn,  ^awi,  that  agree  precisely  with  those  published  by  Mote  for 

the  clamped- free  disk  (although  he  uses  a  different  normaliza¬ 
tion)  (3:  332-333)  .  The  cosine  mode  coefficients  and  are 

determined  by  the  same  process  and  the  second  4X4  matrix.  The  two 
modes,  now  fully  determined,  can  be  evaluated  for  their  dependence  on 
the  torsional  spring  constant  K. 

Consider  the  following  example:  let  a/b  -  0.1,  m  -  0,  and  n  -  2  . 
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The  eigenvalue  ratio  02/fii  this  exaaple  is  plotted  in  Figure  3  as  a 
function  of  the  non* dimensional  parameter  X/D,  the  ratio  of  rotatory 
spring  stiffness  to  disk  bending  stiffness. 


Figure  3 .  Vs .  X/D  for  Histuned  Disk 

a/b  -  0.1  ;  Mode:  m  -  0,  n  -  2 

From  the  figure,  it  is  clear  that  the  eigenvalues  converge  in  the 
limit  of  X  0 ,  and  it  is  clear  from  the  discussion  in  the  previous 
section  that  the  eigenfunction  coefficients  will  converge  as  well,  since 
the  boundary  conditions  for  the  two  modes  become  identical  as  X 
approaches  zero. 

Figure  3  also  illustrates  the  flexibility  of  the  model.  For  this 
example,  the  percent  of  mistuning  (^(wj  —  wi)/wj)  can  be  varied  over  a 
range  of  0  -  62.6%  corresponding  to  the  spring  constant  range  X  -  0  -  *. 
The  maximum  amount  of  mistuning  available  in  the  model  is  far  greater 
than  the  amount  needed  to  model  a  typical  manufactured  disk.  This  is 
fortunate  because  the  mistuned  mode  shapes  will  differ  markedly  for 
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K/D  >  1,  and  the  model’s  ability  to  correctly  model  generic 
imperfections  would  become  questionable. 

Thus,  the  mistuning  model  developed  in  this  section  appears  to  be 
both  a  reasonable  and  flexible  model  of  the  mistuning  phenomenon  in 
general,  although  not  an  accurate  model  of  any  particular  type  of 
imperfection.  To  confirm  the  model's  precision,  the  model's  ability  to 
predict  rotating  node  lines  observed  experimentally  in  stationary 
mistuned  disks  must  be  examined. 
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IV.  stationary  QjLaJs  Forced  Esaana.^ 

The  objective  In  developing  the  expression  for  the  forced  response 
is  to  detemlne  if  it  can  have  rotating  node  lines .  The  non-homogeneous 
differential  equations  of  motion  will  be  developed  and  solved  for  the 
disk  response.  Then  the  conditions  necessary,  if  any,  for  the  response 
to  have  rotating  node  lines  will  be  considered.  A  stationary  point  load 
oscillating  in  time  like  that  used  in  the  experiment  reported  by  Stange 
and  MacBain  will  be  used  first  (1:S).  Then  a  rotating  point  load  will 
be  used  because  it  is  similar  to  a  stationary  point  load  exciting  a 
rotating  disk. 

General  Forced  Response 

According  to  the  expansion  theorem,  as  referenced  in  Chapter  II, 
the  forced  response  can  be  written  as  the  sum  of  the  solutions  to  the 
homogeneous  equation  of  motion  (9:143).  The  expansion  is  expressed  as 


where  the  Wt's  are  the  eigenfunctions  derived  as  solutions  to  the 
homogeneous  equation  of  motion,  and  the  aj(t)'s  are  the  time  dependent 
modal  coefficients. 

The  response  of  a  slightly  mistuned  disk  excited  near  a  resonance 
can  be  approximated  by  the  superposition  of  the  responses  of  the  two 
modes  associated  with  the  two  closely  spaced  resonant  frequencies,  as 
long  as  the  disk  has  no  other  resonant  frequencies  in  this  range.  This 
approximation  is  essential  to  simplifying  the  math  associated  with 
calculating  the  forced  response.  Using  this  approximation,  the  disk 
response  can  be  expressed  as  follows: 
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u;(r,0,O  ■  ai(O/i(r)sinn0  +  a2(O/2('‘)cosn0 


(35) 


where 


and 


Response  A  Stationary  1^2A4 

A  stationary  point  load,  such  as  that  used  In  the  experiments 
reported  by  Stange  and  MacBaln  (6:5),  can  be  represented  as  follows: 


<7 ( /■ .  0 . 0  “  <7o <5( 0  )6(^  - ) Sin  60 f 

where 

qg  -  the  magnitude  of  the  load 

-  the  radial  coordinate  of  the  load 
6  ••  the  Dirac  Delta  function 

w  -  the  frequency  of  the  load 

-  the  angular  coordinate  of  the  load 


(36) 


Undamped  Response.  The  undamped  disk  response  is  determined  as 
the  solution  to  the  non -homogeneous  differential  equation  of  motion. 

The  equation  of  motion  for  the  disk,  given  by  equation  5,  was  derived 
assuming  an  undamped  structure.  Therefore,  the  assumed  response,  given 
by  equation  35,  and  the  load,  given  by  equation  36,  can  be  substituted 
into  equation  5  and  the  time  dependence  solved  for.  Performing  these 
substitutions  yields 


D{a,(OV'‘[/,(r)sinn0]  +  a2(OV^[/2(r)cosn0]} 

+  p/i[d  i(t)/  i(r)sinn9  +  d2(t)/2('')cosn0] 

-qo6(0-0,)6(r-r,]sincot 
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Dateraination  af  Response .  Multiplying  equation  37  by  each 
mode  and  integrating  both  resulting  eqtiations  over  the  domain  of  the 
disk  will  result  in  two  independent  second  order  ordinal^  differential 
equations.  Their  Independence  will  be  due  to  the  orthogonality  of  the 

mistuned  modes.  In  essence,  multiply  by  a,(f)/,(r)sinnfl  and 
ci2(O/2('’)cosn0  and  integrate  to  obtain 


+  ( r,)  sin  sin  to  f 

mi 


cijCO-^  2(0  ”  sin  to/ 

m2 


where  the  modal  masses  and  natural  frequencies  are  given  respectively  by 


m,-‘nph  f  f^ir)rdr 

J  a 


ph 


To  solve  for  the  unknown  modal  coefficients,  a^it)  and,  SjCt),  a 
solution  of  the  form 


a,(/)  -  /I, sin  to/ 


is  assumed  and  then  substituted  into  equations  38  and  39  to  solve  for 
the  amplitudes 


<7o^<;/i(r  Jsinng, 
fc  I  ( 1  -  co^/to?) 
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^  _<?o'’,/a(r,)cosne, 
^  A:a(l -cu^/coj) 

where  the  modal  stiffnesses  are  given  by 

/Cj  -  /rijCof 


(44) 


(45) 


Substituting  these  results  Into  equation  42  and  the  resulting  form  of 
equation  42  Into  equation  35  yields 


w(r.e,t) 


/,(r,)sinnfl, 
/Ci(  1  -cu^/co?) 


/,(r)sina0 


/a(r,)cosree, 
AC2(  1  -  co^/col) 


/2(r)cosn0 


<7or,sinco( 


(46) 


Behavior  Nodes .  The  undamped  response  to  a  stationary 
load  (eq.  46)  does  not  predict  a  rotating  node  line.  To  see  this, 
consider  a  given  radial  coordinate  and  set  the  response  equal  to  zero  as 
shown  below: 


/i(^t]/i(*'o)sinne^ 
/Ci(  1  -ou^/ou?) 


sin  n0 


nod* 


/2('‘,)/2(ro)cosa0, 


cos  n6 


nod9 


1  -cu^/cuj) 

where  rg  -  the  radial  coordinate  of  the  node  point 


qTgr  sincof-O  (47) 


(Note  that  there  will  be  2n  node  points,  n  being  the  number  of  nodal 
diameters  in  the  given  mode.)  Since  the  time  dependence  is  not  zero  for 
all  time,  the  condition  which  must  be  satisfied  is 
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( 1  -co^/co?) 


sinn0„„^,+ 


(48) 


where  X  and  Y  are  constants. 


Solving  for  OiMd.  gives 


®  nod* 


1 

-tan 

n 


(l  -a)Va>?)K 
( 1  -  oj^/cojjA’ 


#  function  of  time 


(49) 


For  a  given  excitation  frequency  u,  the  node  points  will  be  stationary. 
As  the  excitation  frequency  changes,  the  node  points  will  nove  around 
the  disk,  but  they  can  only  have  a  specific  location  for  a  given 
excitation  frequency  u. 


Ih£  Distinction  Between  Node  Lines  and  Node  Points. 

Equation  49  applies  only  to  a  given  radial  coordinate  Tq,  but,  if  the 
mistunlng  is  small,  will  be  approximately  independent  of  r^.  This 

is  because  for  all  Tq  in  a  slightly  mistuned  disk, 


/l('-o)-/2('-o) 


(50) 


Thus,  the  node  points  will  form  a  nearly  straight  node  line  in  such 
cases.  Henceforth,  when  applying  exact  mathematical  equations  to  a 
specified  radial  coordinate  rg  where  a  node  is  desired,  the  term  "node 
point"  will  be  used;  whereas,  when  referring  only  to  the  gross  disk 
behavior,  the  term  "node  line"  will  be  used  instead  with  the 
understanding  that  the  line  is  not  necessarily  straight. 

Damped  Response .  Since  the  undamped  disk  did  not  yield  rotating 
node  lines,  a  slightly  more  realistic  model  is  needed.  The  equation  of 
motion  for  the  damped  disk  cannot  be  readily  derived;  however,  some 
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..*•1;  vS" 


viscous  nodal  daaplng  can  be  easily  added  to  the  independent  second 
order  equations  38  and  39  to  get  an  idea  of  the  possible  effects  of 
danping.  If  the  damping  is  light,  the  eigenfunctions  obtained  from  the 
vmdanped  equation  of  notion  will  not  be  greatly  affected,  since  they 
must  converge  to  the  undamped  form  in  the  limit  of  zero  damping.  As  it 
turns  out,  the  introduction  of  viscous  damping  into  the  equations  of 
motion  does  lead  to  the  existence  of  rotating  node  lines. 

Determination  of  Response .  Recasting  equations  38  and  39 
with  viscous  damping  yields 


d,(0-^2f  ,co,d,(0  +  c«?ai(0-  ^^/i{'*,)sinn0,sina)f 

m  1 

and 

<22(0  +  2^2CU2d2(0-^««2a2(0“  ^^^/j(r,)cosne  sincof 

m2 

where  the  are  the  modal  fractions  of  critical  damping. 


(51) 


(52) 


To  solve  these  equations,  a  solution  of  the  following  form  is 
assumed: 


a,(0  -  /l,sin(co(-a,) 


(53) 


Substituting  this  equation  into  equations  51  and  52  and  solving  gives 


and 


A 


_ <7or<,/i(r  Jsin/ze, _ 

A:,[(l  -ajVaj?)^  +  (2Cia>/a;,)^]" 


tana, 


2^ ,  ou/cu , 
1  -coVoj? 


Substituting  equation  54  back  into  equation  53  yields 


(54) 


(55) 
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qor,/i(r,]sina0,  ,  ^  ^  ^ 

a  ,  (O  - - j  sin  [cut  -  a ,  J 

A:i[(l  -  tt>^/co?]^  +  (2f, co/co, 


Similarly  for  mode  2, 


‘7o'',/2(r,)cosn0, 

1 

A:2[(1  -  (V^  /  ojI]^  +  {2^  2<^  ^  2V]‘ 

2^200/002 

tana, - 2 - - 

l-coVco? 


qo'-,/2('-Jcosne,  .  ^  ^  ^ 

0,(0  - - fSin[a;(-  a,) 

/c2[(l-ou^/a)|)^  +  (2^2t*^/<^2)^]* 


Using  trigonometric  identities,  equations  56  and  59  can  be 


expressed  as 


a,(0 


_ <?o^<7/i(r  Jsinn0, _ 

A:i[(l  -  ou^/cu^]^  +  {2^  icu/uj  i^y 


-(sincotcosa, -cosGofsina,)  (60) 


02(0 


<7o'',/2(r,)cosn0, 


-(sin  00  (cos  a,  "COS  cut  sin  a  2  )  (61) 


A:2[(  I  - cu ^/cu 2)^  +  (2^2f^/cu 2)^  * 


These  results  are  then  substituted  into  equation  35  to  yield  the 
response  given  below: 
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wir 


_ Qorgf^{^^]sinn^^ _ 

A:,[(  1  -oj  Vco?)’  +  (2^,a>/a>,)^]' 


[cosa ,  (sin  cot  sin  nfl)  -  sin  a ,  (cos  cot  sin  n0)]/ ,  (r) 

^ _ qfor,/2(rJcosnfl<, _ 

A:2[(1  -co^/coj)^-*- (2^300/(02)^]* 

[cosa2(sin  cotcosn0)  -  sin  a2(cos<otcosn6)]/2(r)  ^ 

Again,  the  response  will  be  examined  for  the  possibility  of  rotating 
node  lines. 

Behavior  of  Nodes .  The  response  of  equation  62  can  have 
rotating  node  lines  under  certain  conditions.  If  attention  is  again 
restricted  to  a  single  radial  coordinate  r;,,  as  in  the  case  of  the 
undamped  disk,  equation  62  can  be  expressed  as 


w{  r  Q,e  ,t]»  R  ,(co)[cosa,(sinco(sinn0)-sina,(cosco£sin  n0)] 

)[cosa2(sin  co£cosn0)  -  sin  a  j( cos  co  £cosn0)] 

where  and  are  the  modal  response  amplitudes  given  by 


R,(cu) 


<^orJx 

(^<7 

)sin  n0. 

*,[ 

( 1  -  cu’/cu?) 

2^  ,  ou/co  1  ] 

|2  ■ 

1/2 

/?2(a))  = 


qo'",/2 

('■<7 

)cosn0. 

k2[ 

(  1  -  CU^/COj) 

2^2^^/“^  2] 

1^; 

1/2 

/  2  ( ^  0 ) 


Notice  from  equation  63  that  if 


-/?,(co)cosa|  =  R2(^)sina. 
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R,(a))sina,  =  R2i(Ji>')cosa2 


(67) 


then  the  response  can  be  written  as 


w/(r(, ,0,f)  -  ,  (co)[cosa , (sin cxjfsin  n.9 cos u> t cos n9) 

+  sina,(sin6o^cosn0-cosco^sinn0)]  (68) 

or,  using  another  trigonometric  identity, 

uj[rQ,9,t)~  R  ,  (cu)[cosa ,  cos(cof  -  n0)  +  sina,sin(cof-n0)] 

Finally,  this  expression  can  be  reduced  to 


where 


u;(ro,0,f)-/?,(cu)sin(a>(-n0  +  y) 


(70) 


tany  =  cota,  (71) 

Equation  70  is  an  expression  of  the  rotating  node  point  being  sought. 
Observe  that  the  node  point  must  rotate  at  the  angular  speed 


d9  cu  (72) 

dt  n 


which  is  in  the  positive  6  direction.  This  case  will  be  referred  to  as 
the  foxnirard  rotating  node  point. 

If,  instead  of  applying  equations  66  and  67,  the  following 
equations  are  applied. 


R  ,  ( CO  )  cos  a 

,  =  )sin  Oj 

(73) 

R  ,  (co  )sin  a  , 

=  -  )C0S  Oj 

(74) 
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Chen  the  response  can  be  written  as 


ii;(r(,,0,f)-/?,(a))sin(a)f+n6  +  y) 


(75) 


This  case  will  be  referred  to  as  the  backward  rotating  node  point 
because 


(76) 

dt  n 

Eqtiatlon  70  (or  75)  applies  only  to  a  given  radial  coordinate  Tq, 
but,  If  the  mlstunlng  Is  small,  the  conditions  expressed  by  equations  66 
and  67  (or  73  and  74)  will  almost  be  satisfied  for  all  radial 
coordinates.  Again,  as  In  the  case  of  the  stationary  node  lines  of  the 
undamped  disk,  the  lines  will  not  be  exactly  straight,  but  in  this  case 
they  are  rotating.  The  nodes  that  satisfy  equations  66  and  67  (or  73 
and  74)  will  be  referred  to  as  rotating  node  points;  however,  the 
corresponding  gross  response  of  the  disk  will  be  considered  as  having 
rotating  node  lines  even  Chough  they  are  not  always  lines  in  the  strict 
mathematical  sense.  Further  consequences  of  approximate  satisfaction  of 
equations  66  and  67  (or  73  and  74)  are  considered  in  the  "Total  Disk 
Response"  subsection. 

Conditions  Required  isJL  &  Rotating  Node  Point.  The 
equations  that  must  be  satisfied  for  a  rotating  node  point  (eqs  66  &  67 
or  73  &  74)  are  not  very  intuitive.  To  better  understand  the 
requirements  for  rotating  node  points,  the  requirements  should  be 
formulated  in  terms  of  the  relationships  between  the  two  modal  phase 
angles  and  the  two  modal  response  amplitudes.  Equations  66  and  67 
represent  two  homogeneous  equations  in  two  unknowns.  Non- trivial 
solutions  are  possible  only  if  the  determinant  of  coefficients  equals 
zero.  The  matrix  form  of  the  equations  is 
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(77) 


cosa,  sina^ 
sina,  -cosaj 

Taking  Che  determinant,  and  setting  It  equal  Co  zero  yields 


(R,(cv)\  (0\ 

\/?2(a>)/“\o/ 


-  cosa  I  cos  a  2  -  sin  a ,  sin  a  2 


0  -  -cos(a ,  -  ttj) 


or 


a,-a2 


n  . 

—  ±  kn 
2 


(fc-O.l  .2....) 


(78) 

(79) 


The  only  possible  solution  to  this  equation  when  only  two  modes  are  used 
to  describe  the  response  and  (4)2  >  Is 

n  (80) 

fc-O  or  a,-a,-  — 

'  ^  2 


Substituting  this  solution  back  Into  equation  66  gives 


-  R  ,  (cu)cos  a  ,  -  R 2^0^ 


(81) 


or 


(82) 


On  the  other  hand,  If  equations  73  and  74  are  satisfied  with 
W2  >  Wj,  equation  78  still  applies,  but  equation  82  becomes 


Thus  the  apparent  requirements  for  a  rotating  node  point  are  that  the 
response  amplitudes  of  the  two  modes  be  equal  in  magnitude  and  that  the 
phase  angles  differ  by  90  degrees. 
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Distinction  Between  Node  Line  Rotation  Directions.  It  is 
not  clear  from  the  previous  subsection  which  way  the  node  line  will 
rotate  in  a  given  case.  Clearly  both  directions  are  possible,  since  the 
response  amplitudes  and  can  be  positive  or  negative 

depending  on  the  value  of  0^.  Thus,  the  node  line  rotation  can  be 
reversed  by  shifting  by  jr/2n  radians.  By  the  same  logic,  the  case  of 
0)2  <  can  be  evaluated.  In  equation  80,  the  condition  W2  >  was 
assumed  because,  for  the  model  developed  in  Chapter  III  and  the  notation 
established  therein,  W2 ^  with  the  equality  being  satisfied  only  for 
K  —  0.  However,  if  the  modes  are  not  oriented  as  in  the  model  of 
Chapter  III  (in  other  words,  if  the  higher  frequency  mode  is  not  the 
cosine  mode)  then  the  inequality  W2  <  *^1  applies  and  equation  80  becomes 


and  equation  83  becomes  its  partner  equation  in  generating  a  forward 
rotating  node  point.  Yet,  if  equation  83  must  be  satisfied  instead  of 
equation  82,  the  load  angle  must  be  shifted  by  n/2n  to  preserve  the 
forward  rotation  direction.  Therefore,  if  the  modes  are  re -oriented 
such  that  u>2  <  but  the  load  angle  is  not  shifted  by  ir/ln  then  the 
node  line  will  rotate  in  the  opposite  direction.  Thus,  reversing  the 
orientation  of  the  modes  with  respect  to  the  load  position  has  the  same 
effect  on  the  node  line  rotation  direction  as  applying  applying 
equations  73  and  74  instead  of  equations  66  and  67,  since  both  cases 
require  a  shift  of  w/2n  in  6^  to  meet  the  rotating  node  point 
requirements . 

An  intuitive  explanation  may  clarify  the  issue  of  node  line 
direction  of  rotation.  Put  simply,  if  some  set  of  rotating  node  line 
requirements  is  met,  the  node  line  will  rotate  in  the  direction  of  the 
nearest  node  line  of  the  mode  with  the  higher  natural  frequency.  It 
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will  do  so  because  the  mode  with  the  higher  natural  frequency  is  not 
lagging  behind  the  load  as  much  In  time  so  the  rotating  node  line  must 
rotate  in  the  direction  of  its  node  line  (i.e.  in  the  direction  of  the 
line  of  peak  response  of  the  mode  with  the  lower  natural  frequency) . 

Henceforth,  the  inequality  Wj  >  will  be  assumed.  Corresponding¬ 
ly,  the  conditions  of  equations  66  and  67  and  their  corresponding  phase 
angle  and  response  amplitude  requirements  (eqs  80  and  82)  will  be 
considered  to  be  the  requirements  for  a  rotating  node  point  in  the 
stationary  disk  excited  by  a  stationary  load,  as  the  consideration  of 
all  possibilities  would  be  too  confusing. 

Total  Disk  Response .  Thus  far,  only  the  conditions  required 
for  a  rotating  node  point  have  been  developed.  The  response  of  the  rest 
of  the  disk  is  not  obvious  at  this  point;  however,  as  shown  previously 
(eq  50),  if  the  mistuning  is  small,  the  total  disk  response  can 
approximate  rotating  node  lines.  This  can  be  clarified  by  illustrating 
how  the  node  lines  move  through  the  disk  due  to  the  phase  difference 
between  the  two  modal  responses. 

Considering  equations  79,  81,  and  50,  one  can  imagine  the 
responses  of  the  two  individual  modes  as  that  illustrated  in  figure  4 
for  the  mode  n  -  1 ,  m  -  0 .  They  have  approximately  equal  amplitude  but 
as  the  first  reaches  maximum  positive  displacement  the  second  passes 
through  zero  displacement.  Observe  from  eqviations  60  and  61  that 
aj^(C)  -  0  when  wt  -  +  kx  (k  -  0,1,2,  .  .  . )  .  At  such  a  time,  perfect 

node  lines  must  occur  at  all  node  lines  of  mode  1.  At  an  instant  later, 
the  displacement  of  mode  1  begins  to  decrease  and  the  displacement  of 
mode  2  increases.  The  nodes  must  then  retreat  from  the  line  of  peak 
response  of  mode  2  in  the  direction  of  negative  displacement  of  mode  1. 
This  is  how  the  node  lines  move.  Although  only  one  node  point  need  move 
at  the  angular  speed  u/n,  the  other  node  points  must  also  move  at  the 
same  average  speed  and  in  the  same  direction.  Indeed,  if  the  two  mode 
shapes  have  equal  amplitude  at  several  radial  coordinates,  all  node 
points  of  such  coordinates  must  rotate  at  equal  angular  speed. 


Coordinates  having  unequal  response  anplltudes  will  move  faster  or 
slower  than  the  angular  speed  w/n,  depending  on  which  mode  has  the 
greatest  response  at  a  given  Instant  of  time.  Thus,  If  the  rotating 
node  point  requirements  (eqs  76  &  78)  are  met,  the  total  disk  response 
will  consist  of  at  least  one  node  point  rotating  at  the  angular  speed 
w/n  and  all  other  points  rotating  alternately  faster  and  then  slower 
than  this  speed. 

The  speed  variations  of  each  node  point  In  the  disk  cause  the  node 
line  to  be  distorted  whenever  both  coefficients  a^Ct)  and  SjCt)  are 


-  -  j 


C^J 


■ 

■ 
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non-zero.  In  traversing  the  node  line  radially  away  from  the  rotating 
node  point  at  rg,  one  must  veer  away  from  the  nearest  diameter  of  peak 
response  of  the  mode  with  the  greatest  response  amplitude  for  the 
particular  radial  coordinate  (remendser  that  the  modal  response 
amplitudes  are  only  equal  at  rg) .  Observe  that  for  modes  in  which 
n  >  1,  the  shape  of  the  node  line  will  reverse  directions  as  it  moves 
around  the  disk.  This  is  because  the  nearest  diameter  of  peak  response 
of  the  two  modes  exchange  places  with  respect  to  the  rotating  node  line 
as  it  traverses  a  nodal  diameter  of  either  of  the  modes. 

The  insight  provided  by  the  intuitive  look  at  the  mistimed  disk 
response  of  this  subsection  leads  one  to  wonder  what  might  happen  if  the 
rotating  node  point  requirements  (eqs  80  &  82)  cannot  be  exactly  met  for 
any  radial  coordinate  rg  of  the  disk.  Such  a  situation  is  the  subject 
of  the  next  subsection. 

Approximate  Satisfaction  of  Rotating  Node  Point  Require¬ 
ments  .  Thus  far  it  has  been  shown  that  the  response  of  a  mlstuned  disk 
can  have  the  form  of  rotating  node  lines  under  very  stringent 
mathematical  conditions;  however,  these  conditions  are  really  more 
stringent  than  necessary.  It  can  now  be  shown  that  the  requirements  for 
rotating  node  lines  are  actually  much  less  stringent.  This  observation 
would  help  one  to  believe  the  model  is  consistent  with  the  rotating  node 
line  response  reported  by  Stange  and  MacBain,  since  it  is  unlikely  that 
the  requirements  of  both  equation  780  and  equation  82  were  met  in  their 
experiments.  Since  the  mistuning  in  their  experiment  was  small,  it  is 
not  likely  that  equation  80  could  have  been  satisfied,  and  it  would  be 
encouraging  to  know  that  it  didn't  need  to  be  satisfied.  It  should  also 
be  noted  that  the  disk  used  in  their  experiments  was  probably  not 
viscously  damped,  since  this  is  an  additional  reason  for  the  predictions 
made  in  this  thesis  to  differ  from  experimental  results. 

In  the  previous  subsection,  it  was  found  that  points  other  than 
those  at  the  specified  radial  coordinate  rg  can  also  move  around  the 
disk  at  an  angular  speed  similar  to  w/n  leads  to  an  important 
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consideration.  Suppose  the  modal  damping  fractions  are  such  that 
equation  80  cannot  be  satisfied  for  any  radial  coordinate  or  excitation 
frequency.  Such  a  condition  would  not  guarantee  that  rotating  node 
lines  could  not  be  observed  experimentally.  This  Is  clear  from  the 
discussion  of  the  previous  subsection  because  only  a  finite  number  of 
coordinates  can  satisfy  the  requirements  of  equations  80  and  82,  and  yet 
all  node  points  on  the  disk  can  move  at  approximately  the  same  angular 
speed.  Thus,  the  requirements  for  experimental  observation  of  rotating 
node  lines  must  be  considerably  less  stringent  than  those  of  eqtiatlons 
80  and  82 . 

Obviously,  If  equations  80  and  82  are  not  satisfied  there  will  not 
be  a  node  point  rotating  at  exactly  u/n,  but,  for  practical  purposes, 

t 

the  exact  speed  of  the  node  point  Is  not  critical.  A  perfect  node  line 
must  still  occur  when  either  mode  reaches  a  point  of  zero  response,  and 
another  perfect  node  line  must  occur  at  another  location  when  the  other 
mode  reaches  its  time  of  zero  response.  The  only  way  the  node  lines  can 
avoid  moving  is  for  the  two  modes  to  be  exactly  in  phase.  Thus,  it  is 
really  only  the  shape  of  the  rotating  node  line  that  is  affected  by  the 
degree  to  which  equations  80  and  82  are  satisfied  (or  that  non-viscous 
damping  is  present).  Certainly,  as  the  disk  falls  further  away  from 
satisfying  these  equations,  there  will  come  a  point  where  the  node  line 
will  be  too  distorted  to  observe.  The  condition  of  approximate 
satisfaction  of  equations  80  and  82  will  be  addressed  further  in  example 
problem  #4.  At  this  point,  the  final  development  of  the  rotating  node 
point  requirements  must  be  addressed. 

Requirements  211  Modal  Damp inz .  The  response  amplitude  and 
phase  difference  requirements  for  rotating  node  lines  of  equations  80 
and  82  can  be  precisely  met  under  specific  conditions.  These  conditions 
can  be  formulated  in  terms  of  the  modal  fractions  of  critical  damping 
and  the  excitation  frequency  w.  This  will  help  the  designer  or 


experimenter  to  determine  whether  or  not  a  given  disk  can  respond  with  a 
rotating  node  lines  and,  if  it  can,  at  what  excitation  frequency  the 
rotating  node  lines  will  occur. 

Treating  the  damping  values  as  unknowns,  the  requirements  for  a 
rotating  node  point  consist  of  two  eqxiations  in  three  unknowns  (fj^, 
and  w) ;  thus,  one  unknown  must  be  specified.  If,  for  instance,  one 
damping  value  is  specified,  the  remaining  unknowns  can  be  solved  for. 

The  explicit  equations  for  the  remaining  unknowns  are  derived  below. 

Implicit  in  equation  55  are  the  following  equalities: 


sin  a, - ^ 

[[\-(X3'^/(XiJY*{2^i(X3/(XiiYY 

and 


cosaj  - 


1  Vco? 


[( 1  -00^/00?)^  + (2^,0) /O),)^]* 


or 


r/,  2,  2'i2  ,  >>2-15  2Cito/a>j  1-ouVcof 

[  l-o;  /a»J  +  2^,02/02,  I  - - ^ 

L  •'  J  <?ln/Tf  COS  a. 


sin  ttj 


■  I 


(85) 


(86) 


(87) 


The  left  hand  side  of  this  equation  is  found  in  the  denominator  of 
both  response  amplitudes  (eqs  64  &  65);  therefore,  this  equation  can  be 
substituted  into  equation  61  to  yield 


/i(r  J/i(/~o)sinn0 
*:i(  1  -cu^/co?) 


<7  2 

-COS  a, 


/2(r,)/2(ro)cosn0 

^2(2^  2t^/<^2) 


<7  2 

-Sin  a. 


(88) 


From  equation  80,  the  following  equation  is  formed: 


sin^Oj  ”  cos^a , 


(89) 
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Applying  this  to  eqiMitlon  82  yields 


’/l('■,)/l('■o)sin/l0,  /2(r,)/a(ro)cosnd,]  .  3 

- 7 - ? - - i - ; —  Sin  a, 

A:,(l -co’/co?)  aco/coj)  J 

since  sin  02  •*  0,  we  snist  have 

'/i(^J/i(^o)sinne,  ^  /3(/-^)/a(ro)cosrte,~|  ^ 

iti(  1 -co^/co?)  k  2(2^  2<*>  /  7)  J 

Solving  this  quadratic  equation  gives 


(90) 


(91) 


/a(''J/2('‘o)/n,coj 


± 


I 

/i('‘J/i(^o)?2*3tanae,Y  ^  ^3 

/2(r,)/2(ro)^ia>2  J  ‘ 


(92) 


Thus,  the  excitation  frequency  necessary  for  a  forward  rotating 
node  point,  can  be  solved  for  if  Is  known.  On  the  other  hand.  If  Ti 
is  known. 


(JO 


frn 


/ 2!^  g)/ 2(^  0)?  I  ^  I 

/i(r,)/i(ro)tann0,m2co 


f  ijl'  q)/2(^  o)Cl^l _ 

,/i(/',)/i(ro)tann0,m2a», 


co: 


(93) 


In  equations  92  and  93 ,  the  positive  root  must  be  taken  to  obtain  the 
desired  positive  frequency. 


36 


Aft«r  on«  aodal  fraction  of  critical  damping  is  asstimed  and  the 
raquirad  excitation  frequency  solved  for,  the  other  modal  damping  value 
can  be  determined  by  rearranging  equation  87.  This  yields 


if  Cj  assumed  and 


(94) 


(95) 


if  Cl  Is  assvimed. 

If  the  modal  damping  of  both  modes  are  assumed  to  be  equal, 
equations  92  and  93  can  be  set  equal  and  the  damping  solved  for.  The 
resulting  expression  is 


(96) 


The  excitation  frequency  corresponding  to  this  fraction  of  critical 
damping  can  be  calculated  using  either  equation  92  or  equation  93  . 

Damping  Limits .  The  amount  of  modal  damping  cannot  exceed 
certain  limits  if  the  rotating  node  point  requirements  (eqs  80  &  82)  are 
to  be  satisfied.  This  limitation  is  due  to  the  phase  requirement  of 
equation  80.  If  the  modal  damping  is  too  high  relative  to  the  amount  of 
mistuning,  the  required  phase  difference  cannot  be  achieved  for  any 
excitation  frequency.  The  maximum  values  that  the  fractions  of  critical 
damping  can  take  on  can  be  calculated,  even  though  the  equal  damping  as 
calculated  in  equation  96  may  be  more  useful. 
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The  precise  llnltetlons  on  the  aodel  fractions  of  critical  dasplng 
can  be  determined  by  considering  that  the  phase  angle  of  either  of  the 
two  modes  is  90  degrees  at  the  natural  frequency.  That  is, 

-•  it/2  and  ■■  s/2  .  Thus,  a  rotating  node  point  can  only 

occur  at  excitation  frequencies  between  the  two  natural  frequencies. 
Therefore,  the  maximum  fractions  of  critical  damping  can  be  determined 
by  setting  equations  92  and  93  equal  to  the  corresponding  natural 
frequencies.  This  results  in  the  following  equations: 


2m,/j(r,)/j(ro)  Vto?  ) 


m,/a(rJ/a(ro)  ( 

"  2mj/,(r,)/,(ro)tann0,V  ^\j 


(97) 


(98) 


If  either  of  the  fractions  of  critical  damping  exceeds  these  values,  the 
other  fraction  cannot  be  ouide  small  enough  to  satisfy  the  phase 
requirement  of  equation  80.  Notice,  though,  that  if  one  fraction  takes 
on  its  maximum  value  the  other  fraction  must  approach  zero.  Since  one 
would  normally  expect  the  two  modal  damping  values  to  be  about  equal, 
the  maximum  damping  values  given  by  equations  97  and  98  aay  not  be  very 
useful. 

Example  Problem  »1.  An  example  problem  will  help  to 
Illustrate  the  equations  developed  thus  far.  Consider  the  following 
parameters  for  a  steel  disk: 
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b  - 

0.13  meters 

X  - 

24.0  N-meters/Radian 

a  - 

0.013  meters 

0.05  meters 

0.3 

a/8  radians 

h  - 

0.003175  ffliters 

n  - 

2 

E  - 

206.844  MPa 

a  - 

0 

fi  - 

7800  kg/cu  meter 

^0  - 

0.05  meters 

Applying  the  methods  described  In  Chapter  III  yields  the  following  data: 


^02  1 

-  18.24018461 

^023 

-  18.32274608 

^02 

-  1 .0 

f  02 

-  1 .0 

Bq2 

-0.02630958 

^02 

-0.02640427 

^02 

-0.20438468 

C02 

-0.19803424 

-0.01698596 

H  oa 

-0.01704926 

-2.68107X10® 

-2.72086X10® 

m, 

-0.09894148 

fTX  2 

-0.09861232 

cw, 

-  1646.13  rad/s 

coa 

-  1661.07  rad/s 

The  radial  dependences  of  the  mode  shapes  are  plotted  in  figure  5.  It 
is  clear  from  this  figure  and  from  the  split  natural  frequencies  that 
the  mistuning  is  small.  This  Justifies  the  ability  of  the  mistuning 
model  to  model  small  imperfections  without  drastically  affecting  the 
mode  shapes. 
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Figure  S.  Exeaple  Problem  #1,  Mlstuned  Redial  Mode  Shapes 
K  -  24.0;  a/b  -  0.1;  Mode:  a  -  0,  n  -  2 

Next,  the  condition  of  equal  modal  damping  Is  calculated  using 
equations  96  and  92.  The  results  are 


004516  cu,„-  1653.59 

Using  equations  97  and  98,  the  damping  limits  are  calculated. 

They  are 

-0.009037  -0.009028 

If  either  of  the  fractions  of  critical  damping  are  higher  than  these 
values,  a  node  point  rotating  at  constant  angular  speed  cannot  occur  at 
the  specified  radial  coordinate  rg  regardless  of  how  small  the  other 
damping  value  becomes.  To  demonstrate  this,  figure  6  plots  the  phase 
angles  for  the  damping  fractions  “  0-®^  • 


■'rt  I  I  t  I  >  r  t  I  I  I  '  '  '  '  '  '  '  '  I  '  '  '  '  '  ^  '  I  '  '  '  '  '  '  '  '  M  '  '  '  '  '  '  ' 

1630  1640  1650  1660  1670 

CJ  (RAOIANS/SEC) 

Figure  6.  Exaaple  Problem  #1,  Modal  Phase  Angles 

fi  -  0.01,  fj  -  0.01 

Equation  80,  the  phase  requirement  for  a  rotating  node  point.  Is  clearly 
never  satisfied  for  the  damping  value  specified  here;  therefore,  a  node 
point  rotating  at  the  constant  frequency  w/n  cannot  occur  for  the 
specified  radial  coordinate  r,. 

Next,  one  fraction  of  critical  damping  Is  specified  and  the  other 
parameters  are  calculated  using  equations  93  and  95.  Letting  0.00400 
yields 

^2  =  0-  00503  /r««  “  1 654 , 44 

The  modal  response  amplitudes  for  these  parameters  are  plotted  In  Figure 
7,  and  the  phase  angles  are  plotted  In  Figure  8. 
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Figure  8.  Example  Problem  #1,  Modal  Phase  Angles 

-  0.00400,  r2  "  0-00503,  Rotating  Node  Point 


Th«s«  flg\ir«s  show  thac  the  conditions  required  for  s  rotating  node 
point  (eqs  80  &  82)  are  indeed  satisfied  at  the  indicated  frequency 
(1654.44  rad/sec). 

The  shape  of  the  rotating  node  lines  can  be  plotted  by  varying  the 
tisie  and  radius  in  eqtiation  62.  If  this  is  done  using  the  aodal  damping 
values  and  excitation  frequency  calculated  earlier,  the  node  line  motion 
shown  in  figure  9  results.  In  this  figure,  the  node  line  distortion  is 
greatly  exaggerated.  Notice  that  when  the  node  line'*  cross  the  natural 
node  lines  of  each  mode,  they  become  straight.  At  such  times,  the 
displacement  of  the  other  mode  is  zero.  The  times  c  that  this  will 
occur  can  be  determined  by  setting  t  equal  to  the  angle  of  a  desired 
modal  nodal  diameter  and  setting  the  response  of  equation  70  equal  to 
zero  (this  will  place  the  node  point  at  the  modal  nodal  diameter  and  the 
time  can  be  determined  as  C  -  (jni  -  7)/a>. 


Figure  9.  Example  Problem  »1,  Rotating  Node  Lines 
Mode:  m  -  0,  n  -  2 


Figure  10  is  an  accurate  plot  of  the  node  line  shape  at  its  maximum 
curvature  and  magnified  100  times  along  the  Y-axis.  If  it  was  plotted 
to  scale  or  even  magnified  10  times,  it  would  appear  as  a  straight  line. 
Note  that  the  large  curvature  at  the  disk  outer  edge  is  expected  because 
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the  mode  2  node  shape  Is  bent  downward  at  the  edge  by  the  torsional 
springs.  In  essence,  the  aaount  of  curvature  Is  a  function  of  the 
difference  between  the  mode  shapes  along  the  radial  direction. 


Figure  10.  Example  Problem  #1.  Node  Line  Shape  Magnified  100  Times 
Mode:  m-0,  n-2  ;  JC-24.0 

Example  Problem  #2-  'Fhe  node  line  shape  can  be  plotted  for  other 
modes  as  well.  By  choosing  a  different  mode,  the  effects  of  mlstvinlng 
on  different  modes  can  also  be  explored.  Consider  using  all  of  the  same 
parameters  used  In  example  #1  except  for  the  number  of  nodal  circles  In 
the  mode.  Instead  of  m  -  0,  let  m  -  3.  The  resulting  radial  mode 
shapes  are  plotted  In  figure  11.  The  mlstunlng  Is  so  slight  that  the 
mode  shapes  are  virtually  Identical.  The  modal  data  for  these  modes  are 
as  follows: 


^23  1 

-99.27284419 

^23  2 

-99.29857727 

^02 

-  1  .0 

^02 

-  1  .0 

5  02 

-0.56873390 

F  07 

-0.56903183 

C  02 

- -.00000492 

^02 

- -.00000489 

^02 

-0.54019380 

^02 

-0.54058599 

-9.39254X10^ 

^2 

-9.40461X10^ 

m, 

-0.03950465 

m,  2 

-0.03951444 

02, 

-  48760.427  rad/s 

0J2 

-48785.709  rad/s 
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Figure  11.  Exaaple  Problem  #2,  Mistuned  Radial  Mode  Shapes 
a/b  -  0.1  ;  Mode:  a  -  3,  n  •-  2 

Note  Chat  the  percent  of  alstunlng  <»(w2  ~  wi)/wi)  Is  much  smaller  in 
this  mode  for  the  same  torsional  spring  constant.  This  makes  good 
intuitive  sense  because  one  would  expect  the  boundary  condition  to  have 
less  influence  on  Che  disk  as  a  whole  for  modes  that  have  more  nodal 
diameters  or  more  nodal  circles.  The  boundary  imperfection  becomes  more 
and  more  isolated  from  most  of  Che  disk  by  the  additional  modal  node 
forms,  and  thus  zone  of  Che  disk  being  influenced  by  the  imperfection 
becomes  smaller.  The  main  effect  of  reduced  mistuning  is  that  the  modal 
fractions  of  critical  damping  must  be  smaller  if  equation  80  is  to  be 
satisfied.  This  can  be  easily  shown  by  assuming  one  damping  value  and 
calculating  what  Che  other  must  be. 
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To  plot  the  aodel  phase  angles  and  node  line  shape,  one  fraction 
of  critical  danplng  aust  again  be  specified.  If  f  -  0.00025,  then 
equations  93  and  95  yield 

^2-0.000268  48773.35 

The  modal  phase  angles  are  plotted  In  figure  12,  and  the  rotating  node 
line  Is  plotted  In  figure  13  during  a  time  of  maximum  curvature.  The 
Y*axls  magnification  In  this  figure  Is  10. 


Figure  12 .  Example  Problem  #2 ,  Modal  Phase  Angles 

-  0.000250,  {-2  "  0.000268,  Rotating  Node  Point 


0.05  0.07 

r  (METERS) 


Figure  13.  Exaiiple  Problem  #2,  Node  Line  Shape  Magnified  10  Times 
Mode:  a-3,  n-2;  JC-24.0 


Observe  that  the  distortions  of  the  node  line  occur  at  the  nodal  circles 
of  the  mode.  These  distortions  may  be  attributed  to  numerical 
Inaccuracy  in  the  region  of  a  nodal  circle.  The  angular  position  of  the 
node  line  cannot  be  accurately  determined  in  a  nodal  region  of  the  mode 
shape.  These  distortions,  in  any  case,  would  not  be  apparent  to  the 
experimental  observer,  since  they  are  near  the  nodal  circles  of  the 
mode.  Because  of  the  suspected  numerical  problem  with  this  node  line, 
another  example  will  be  considered  without  nodal  circles. 

Example  Problem  #2.  Now  consider  using  all  of  the  same 
parameters  used  in  example  #1  except  for  the  number  of  nodal  diameters 
in  the  mode.  Instead  of  n  —  2,  let  n  —  4.  The  resulting  radial  mode 
shapes  are  plotted  in  figure  14. 
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Figure  14.  Exeople  Problen  #3,  Mlstuned  Radial  Mode  Shapes 
a/b  -  0.1;  Mode:  a  -  0,  n  -  4 


The  aodal  data  for  these  aodes  are  as  follows: 


^40,  -35.94503247 

^02  —  1  •  0 
Bo2  "  .000021536 
Coa  -0.04279918 
Doa  -  0.00001409 
k^  -0.19118X10' 
m,  -0.04678186 
CO,  -6392.700  rad/s 


^403  -35.98051624 

£o2  -1-0 

Foa  -  .000021649 
Coa  -0.04216615 
//o2  -0.00001416 
itj  -0.19179X10' 
mj  -0.04674543 
coj  -6405.328  rad/s 


Letting  -  0.00009  yields 

^j-0. 001074  a>,,„-6399,54 

The  shape  of  the  rotating  node  line  is  plotted  in  figure  15  magnified 
100  times. 


Figure  IS.  Example  Problem  #3,  Node  Line  Shape  Magnified  100  Times 
Mode;  m-0,  n-4  ;  JC-24.0 

This  node  line  looks  more  like  it  is  expected  to.  The  curvature  is  less 
than  that  of  example  #1,  because  the  percent  mistiming  is  smaller.  The 
mlstuning  is  smaller  because  of  the  greater  number  of  nodal  diameters  in 
the  mistuned  modes.  Again,  the  greatest  curvature  occurs  where  the  two 
mode  shapes  differ  the  most  -  at  the  outer  boundary.  One  more  node  line 
will  be  plotted  later,  in  example  #8,  as  a  more  important  idea  must  now 
be  explored. 


RTc«iiip]ft  Problem  #4.  As  an  example  of  what  might  happen  if 
equations  80  and  82  are  not  satisfied,  consider  the  following  changes  to 
example  #1;  let  Tj  “  0.015  and  fj  ”  0*017  .  These  values  are 
considerably  greater  than  the  maximum  values  allowed  if  a  node  point 
rotating  at  the  constant  frequency  w/n  is  desired;  however,  the  previous 
discussion  about  approximate  satisfaction  of  the  rotating  node  point 
equations  suggests  that  a  rotating  node  line  will  still  occur  in  some 
shape  or  form.  The  phase  requirement  cannot  be  met  with  these  damping 


v«lu«s  as  Is  claar  from  figure  S;  however,  the  nodal  responses  aust  be 
equal  for  soae  excitation  frequency.  Solving  equation  82  for  this 
frequency  and  calculating  the  corresponding  phase  angles  yields 

CO  -  1 659 . 48 1  rad/sec  a ,  -  2 . 06459  radians  a  j  -  1  . 5 1 449  radians 

Note  that  the  phase  requirement  (eq.  80)  Is  far  fron  being  satisfied. 

If  the  frequency  calculated  above  Is  substituted  into  equation  62 , 
the  disk  response  can  be  plotted  for  any  tine,  radius,  or  angle.  By 
plotting  the  total  disk  response  as  a  function  of  r  and  i ,  the  shape  of 
the  disk  can  be  observed.  In  this  way,  one  can  compare  more  than  just 
the  position  of  the  node  lines  -  one  can  compare  the  entire  disk  shape. 
This  will  assist  the  experimenter  In  comparing  his  observations  to  the 
model  predictions.  The  response  Is  plotted  in  figures  16-18  for  three 
<ilffcrent  times  with  radius  as  a  parameter. 


r-.13 

/TioN 
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Figure  16.  Example  Problem  #4,  Disk  Response  Vs.  Theta 

C  -  -0.980  milliseconds;  Ti  “  0.015,  "  0-017 
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Figure  17.  Example  Problem  #4,  Disk  Response  Vs.  Theta 

C  -  -0.897  milliseconds;  -  0.015,  "  0.017 


Figure  18.  Example  Problem  #4,  Disk  Response  Vs.  Theta 

t  -  0.522  milliseconds;  -  0.015,  ”  0.017 
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These  figures  show  thec  the  node  line  is  virtxielly  straight  and 
that  the  total  disk  response  has  the  saae  shape  as  that  of  a  single 
■ode.  Thus,  even  when  only  equation  82  is  satisfied,  the  disk  has  a 
rotating  node  lines  and  has  the  same  general  shape  as  that  of  the  tuned 
mode  shape.  Observe  however  that  the  amplitude  of  the  displacement 
changes  as  the  node  line  moves  around  the  disk.  Since  equation  80  is 
not  satisfied,  the  modes  need  not  add  up  to  the  same  amplitude  for  all 
times.  In  addition,  the  node  line  cannot  rotate  at  the  constant  speed 
w/n  because  there  is  not  enough  phase  shift  between  the  response  of  the 
two  modes.  The  actual  node  line  shape  is  plotted  in  figure  19.  As  one 
might  expect  because  of  the  percent  of  mistuning  (the  same  as  for 
example  #1),  it  looks  much  like  the  node  line  of  example  #1. 


Figure  19.  Example  Problem  #4,  Node  Line  Shape  Magnified  100  Times 
Mode:  m-O,  n-2  ;  Jk-24.0 

The  position  of  the  node  line  is  plotted  versus  time  in  figure  20. 
Note  that  it  moves  faster  and  then  slower  than  the  speed  w/n  but  must 
have  an  average  angular  speed  equal  to  w/n  because  each  mode  must  pass 
through  zero  response  twice  per  cycle  of  excitation. 
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Figure  20.  Exanple  Problem  #4,  Node  Line  Poslclon  Vs.  Time 
-  0.015,  f2  “  0.017;  Mode:  a  -  0,  n  -  2 

A  more  direct  observation  about  the  amplitudes  In  figures  17  and 
18  can  now  be  made.  Observe  that  the  response  of  figure  17  occurs  at  a 
time  when  the  node  line  Is  moving  fast.  This  Is  a  time  when  the  energy 
of  the  disk  Is  mostly  kinetic,  and,  therefore,  the  strain  energy 
(deformation)  of  both  modes  Is  small.  In  contrast,  the  response  of 
figure  18  occurs  when  the  disk  energy  Is  mostly  potential,  and, 
therefore  the  response  of  both  modes  Is  large.  This  explains  why  the 
total  response  In  figure  17  Is  smaller  than  that  In  figure  18,  as  the 
sum  of  two  small  modal  responses  must  be  smaller  than  the  sum  of  two 
large  modal  responses.  Thus,  when  the  node  line  rotates  slowly,  the 
response  Is  large  and  vice-versa.  The  response  In  figure  16  occurs  when 
the  response  of  mode  2  Is  zero;  thus.  Its  amplitude  represents  the  mean. 


The  stationary  alstuned  disk  response  to  a  stationary  oscillating 
load  has  now  been  extensively  analyzed.  The  response  was  determined  as 
were  the  conditions  required  for  rotating  node  lines  to  occur.  The 
conditions  required  for  a  node  point  rotating  at  the  constant  angular 
«peed  w/n  are  that  the  modal  response  astplitudes  be  equal  in  magnitude, 
that  the  modal  responses  be  shifted  in  phase  by  ir/2  radians,  and, 
obviously,  that  the  modes  be  damped  to  cause  this  phase  shift.  It  was 
also  shown  that  as  long  as  there  is  some  phase  difference  between  the 
modal  responses  and  the  modal  response  amplitudes  are  equal  the  node 
lines  will  rotate  albeit  not  at  the  constant  speed  u/n.  Most 
importantly,  it  was  shown  that  analytic  models  can  predict  the  existence 
of  rotating  node  lines  observed  experimentally  by  Stange  and  MacBain. 

The  stationary  load  has  provided  excellent  insight  into  the 
response  of  a  mis tuned  disk;  however,  most  aerospace  applications 
Involve  rotating  disks.  The  next  step  to  be  taken  is  to  simulate  a 
rotating  disk  by  rotating  a  load  around  a  stationary  disk. 

R^SPgPge  A  £2  Load 

A  point  load  rotating  around  a  stationary  disk  at  constant  angular 
speed  can  be  used  as  a  model  for  a  rotating  disk  excited  by  a  stationary 
point  load.  The  resonant  frequencies  of  the  rotating  disk  are  known  to 
be  higher  than  those  of  the  stationary  disk  but  the  mode  shapes  are 
similar.  The  response  of  the  stationary  disk  to  a  rotating  load  will  be 
determined  and  later  considered  as  a  model  of  the  rotating  disk 
response.  The  conditions  required  to  form  rotating  node  lines  will  also 
be  determined.  These  conditions  should  be  slightly  different  from  those 
of  the  stationary  load  case  because  the  two  mistimed  modes  are 
automatically  excited  out  of  phase  by  the  rotating  load. 

A  rotating  load  of  constant  amplitude,  angular  velocity,  and 
radial  position  can  be  described  by  the  following  expression; 


q(r  ,9,  t)~  q  06(9-0}  t)6{r-r^] 


(99) 
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wh«r« 


qg  -  the  oagnltude  of  the  load 

-  the  radial  coordinate  of  the  load 
S  -  the  Dirac  Delta  Function 
w  -  the  angular  speed  of  the  load 

TTnrfai^^d  Response .  Incorporation  of  the  rotating  load  Into  the 
eqviatlons  of  motion  for  the  stationary  load  only  affects  the  force  side. 
Equations  38  and  39  become 


cI,(0-^co?a,(0 


^^/.(rjsinncot 


(100) 


and 


tn  2 


(101) 


A  solution  of  the  form 


a,(0"  /^iSincuf 


(102) 


Is  again  assumed  for  mode  1;  however,  a  different  solution  Is  assumed 
for  mode  2.  It  Is 


ajCO  “  /1 2  cos  CO  I 


(103) 


Different  solution  forms  for  the  two  modes  are  assumed  because  the  two 
mlstuned  modes  are  excited  out  of  phase  by  the  rotating  load.  This 
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should  lead  to  a  change  in  the  conditions  required  to  excite  a  rotating 
node  point,  because  the  phase  of  the  response  is  critical  to  exciting 
such  a  response  (recall  eq  80). 

Determination  of  Response .  Solving  the  independent  differ¬ 
ential  equations  leads  to  the  response,  given  by 


w{r  ,e,t) 


Jfc,(  1  -n^aj^/co?) 


/,(r)sinn0sin  ncot 


Qor^fiir^] 


)t2(  1 


jxf 2ir)cos nd cos n(vt  (104) 


Equation  104  indicates  that  the  disk  vibrates  at  the  frequency  nu, 
which  equals  the  excitation  frequency  ut  only  for  modes  with  n  -  1  (one 
nodal  diameter) .  This  unusual  response  is  the  result  of  using  a  moving 
load.  Such  a  load  will  pass  2n  modal  nodes  in  one  trip  around  the  disk 
and  will  thus  require  only  1/n  of  the  frequency  of  the  stationary  load 
to  excite  the  same  mode.  In  plots  of  disk  response  ,  where  this 
peculiarity  may  cause  confusion,  the  excitation  frequency  and  response 
frequency  will  be  distinguished  by  subscripts. 

Behavior  of  Node  Points .  Again  a  radial  coordinate  must  be 
specified  to  solve  for  a  rotating  node  point.  If  this  is  done,  the 
response  becomes 


qor«,/i(rJ/,(ro) 
A:i(  1  -n^co^/oj?) 


sin/i0sinrtcut 


1  - 


cosndcosncot 


(105) 


which  requires  only  one  condition  to  form  a  rotating  node  point,  and 
that  condition  is  as  follows: 
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(106) 


qor^fi{r^)ft[ro)  <7or,/a(r,)/a(ro) 

ilc,(  1 -n^to^/co?)  Jtj(  1 -n’co^/coi) 

If  this  condition  is  satisfied,  the  response  can  be  written  as 


(107) 

(108) 

which  is  a  node  point  rotating  in  the  positive  i  direction.  On  the 
other  hand,  if 


u;(r<,,0,t)  -  '^^°j(sinncofsin/zg->-cosncotcos/zO) 

1 -n^oj^/cotJ 


or 


u/(ro.0.0’ 


qo/-,/i(r,)/,(ro) 

it,(  1  -n’co*/a>  i) 


cosn(_(vt  -  6) 


qorJi{r^]f^{ro] 
A:i(  1 -n^co^/co?)  fc2(  1 


the  response  would  be 


(109) 


it'(ro.0.O 


qo/’,/i(rJ/,(ro) 

A:,(l  -n’ojVcu?) 


cosn(co(  +  0) 


(110) 


which  is  a  node  point  rotating  in  the  negative  9  direction.  Thus,  in 
the  undamped  disk  excited  by  a  rotating  load,  both  forward  and  backward 
rotating  node  lines  are  theoretically  possible. 

The  rotating  node  line  angular  speed  is  slightly  different  for  the 
rotating  load  case.  Equations  108  and  110  indicate  that  the  node  line 
angular  speed  is  given  by 


dt 


±UJ 


(111) 
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Thus,  the  node  line  angular  speed  due  to  a  rotating  load  Is  independent 
of  the  number  of  nodal  diameters  in  the  mode.  In  fact,  in  the  case  of 
the  forward  rotating  node  point,  the  node  lines  would  rotate  right  along 
with  the  load  (dt/dt  -  +w) . 

Equation  109  implies  that  the  response  of  the  two  modes  is  equal 
in  amplitude  but  opposite  in  sign,  rather  than  equal  in  amplitude  and 
equal  in  sign,  as  in  equation  106.  Equation  109  will  normally  be 
satisfied  somewhere  in  the  middle  of  the  split  natural  frequencies. 
Equation  109  can  be  solved  for  this  excitation  frequency  required  to 
excite  a  backward  rotating  node.  The  result 


J/i(ro)/m,-a>^/2(r  J/2(ro)/ma 


Equation  106 ,  on  the  other  hand,  will  normally  only  be  satisfied 
for  excitation  beyond  the  dual  resonant  frequency  band,  if  at  all. 
Solving  equation  106  for  the  excitation  frequency  needed  to  excite  a 
forward  rotating  node  gives 


<^1/ ijr,)/ i{ro]/m,->-u)U2{r^]/2{ro)/m2 
n^(/i(r,)/,(ro)/m,  +/2(r,]/2(ro)/m2)  . 


The  frequencies  calculated  using  equation  113  may  easily  be  too 
high  to  obtain  an  accurate  solution  because  the  assumption  of  equation 
35  (that  the  disk  response  can  be  approximated  as  the  superposition  of 
the  responses  of  the  two  mistuned  modes)  is  invalid  except  at 
frequencies  near  the  natural  frequencies  of  the  two  mistuned  modes.  For 
this  reason,  a  rotating  node  point  would  not  be  expected  at  this 
frequency,  and  equation  113  should  not  be  used  to  predict  a  rotating 
node  unless  the  modes  of  interest  are  known  to  dominate  the  response  at 
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Che  indicated  frequency.  There  is  no  guarantee  that  eqxiation  113  will 
even  yield  a  real  frequency,  since  the  response  aoplitudes  do  not  have 
to  cross-over  at  any  frequency  above  the  higher  natural  frequency.  The 
amplitudes  can  cross-over  if  Che  modal  parameters  and  load  coordinate  r, 
are  right,  as  can  be  shown  by  altering  Che  mistiming  model  or  moving  Che 
load.  This  will  be  demonstrated  in  example  problems  5  and  6. 

Example  Problem  #5 .  An  example  will  help  illustrate  the 
conditions  needed  Co  generate  rotating  node  lines  in  Che  case  of  Che 
imdamped  disk  with  a  rotating  load.  Consider  using  Che  same  parameters 
as  in  example  #1  except  for  the  coordinate  of  Che  load.  Let  r,  -  0.13, 
Che  outer  radius  (edge)  instead  of  0.05.  Applying  equation  112  yields 
Che  required  excitation  frequency 

CO  -  826  .805  radians/sec 

This  frequency  corresponds  Co  the  node  point  rotating  backward  in  the 
disk  at  the  constant  angular  speed  w  (according  to  equation  111). 
Equation  113,  on  the  other  hand,  does  not  yield  a  real  frequency.  Thus, 
for  the  stated  parameters,  the  modal  response  amplitudes  will  never  be 
equal  with  the  same  phase.  Figure  21  plots  the  response  amplitudes  of 
Che  two  modes  in  a  frequency  band  where  they  might  be  expected  to  cross. 

Although  the  modal  response  amplitudes  are  apparently  never  equal 
at  this  radial  coordinate ,  they  may  be  equal  at  another  coordinate  for 
some  frequency.  Review  of  equation  113  suggests  that  some  alteration  of 
the  modal  masses,  the  coordinate  of  the  load,  or  the  coordinate  of  the 
node  point  may  lead  to  a  real  frequency:  however,  only  the  third  option 
does  not  involve  a  change  in  the  configuration. 
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Figure  21.  Example  Problem  #5,  Response  Amplitudes 

No  Damping,  Tq  -  0.05,  No  Forward  Rotating  Node 


Example  Problem  #6.  Consider  the  same  configuration  as  in 
example  #5,  but  now  let  rg  -  0.125.  Application  of  equation  112  yields 
the  excitation  frequency  required  for  a  bacicward  rotating  node 
cu  =  826  .815  radians/sec 

This  frequency  is  very  close  to  that  of  example  #5,  which  is  simply  an 
indication  that  the  mistuning  is  small.  Equation  113,  however,  now 
gives  a  real  excitation  frequency.  It  is 

CO  =  2 . 093 .415  radians/sec 


Thl*  is  the  frequency  required  for  the  forwerd  rotating  node  point  in 
the  undaiiped  disk  subjected  to  a  rotating  load.  Thus,  the  point 
r,  -  0.05  cannot  rotate  at  the  excitation  frequency  but  the  point 
r,  -  0.125  nay.  Figures  22  and  23  plot  the  nodal  response  amplitudes 
for  the  paraneters  of  this  exanple  in  two  different  frequency  bands  • 
where  eqxiations  106  and  109  respectively  are  satisfied.  The  difference 
between  the  nodal  response  amplitudes  is  also  plotted  in  figure  23 
because  the  two  amplitudes  are  too  close  to  observe  their  crossing. 
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Figure  22.  Example  Problem  #6,  Response  Amplitudes 

No  Damping,  Tg  -  0.125,  Backward  Rotating  Node 

These  figures  clearly  show  that  the  response  amplitudes  of  the  two 
mlstuned  modes  at  the  frequency  required  for  the  forward  rotating  node 
point  is  too  small  for  these  modes  to  be  considered  dominant. 

Therefore,  a  forward  rotating  node  point  would  not  really  be  expected  at 
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the  Indicated  excitation  frequency  of  2,093.415  radians/sec.  To  be 
certain  of  this,  the  other  resonant  frequencies  of  the  disk  would  have 
to  be  determined. 


Wr..p<ms*  (RADIANS/SEC) 


Figure  23.  Example  Problem  #6,  Response  Amplitudes 

No  Damping,  rg  -  0.125,  Forward  Rotating  Node  (?) 

The  undamped  disk  has  provided  some  interesting  results  in  the 
case  of  the  rotating  load;  however,  there  is  no  such  thing  as  an 
undamped  disk.  The  effects  of  damping  on  the  previous  results  must  be 
addressed. 

nampod  Response.  To  determine  the  damped  response  due  to  a 
rotating  load,  the  damping  is  again  added  to  the  undamped  disk  equations 
of  motion.  The  resulting  equations  are  solved  for  the  modal  responses, 
and  those  responses  are  added  by  employing  the  expansion  theorem.  The 
response  will  again  be  examined  for  the  possibility  of  rotating  node 


lines  and  the  conditions  needed  to  achieve  then. 

If  viscous  nodal  danping  is  again  added  to  the  differential 
equations  governing  the  response,  a  change  in  the  rotating  node  point 
requirements  would  be  expected.  This  is  in  fact  what  happens,  as  will 
be  sho«m. 

Adding  viscous  danping  to  equations  100  and  101  yields 


mi 


and 


2COjdj(0-*‘C0  2aj(0“^^/2('',)cosrtcoi 

m.2 


(114) 


(115) 


For  mode  1,  a  solution  of  the  fora 

a,(()-  /4,  sin  (/ICO  a,) 

is  assumed;  whereas  for  mode  2,  the  assumed  solution  is 


(116) 


<^2(()  ■ '^acos(/ico(-aj)  (117) 

Determination  of  Response.  Substituting  equations  116  and 
117  into  the  differential  equations  of  notion  (114  and  115)  and  solving 
for  the  response  at  a  specified  radial  coordinate  gives 

u;(ro  ,0 ,()  -  R ,  (co  )[cosa ,  (sin/icotsin  /i0)-  sin  a ,  (cos /i  cot  sin  n0)] 

+  )[cosa2(cos/icotcos/i0)  +  sina2(sin/icot cos/i0)]  (118) 

where  i?i(c>))  and  are  again  the  modal  response  amplitudes  given  now 

by 
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Qor^fi[r,]fi[ro] 


fci[(  1  -n^co’/co  i)*  +  (2f  ,n( 
<7or,/2(r,)/2(ro) 


/C2[(i  -  n^cv^/ojiy 


The  phase  angles  are  now  given  by 


2n^jaj/ajj 

tana.  - ; — r - r 

1  -  n^oj^/ajf 


(121) 


Now,  if 


/?  ,(a))cosa ,  -  /?  2(00)  cos  a  2  =  0 


/?  i(a))sin a,  - /?2(<«)sina2  =*  0 


(122) 


(123) 


equation  118  can  be  written  as 


/(rg,0,f)  =  R^  (  CO  )  cos  a  ,cosa(cof-0)  +  ^,(co)sina|Sinn(co(-0)  (124) 


where 


tan  y  =  cota , 


(125) 


(126) 


Equation  125  is  an  expression  of  the  forward  rotating  node  point. 

Notice  that  the  bacicward  rotating  node  point  Is  no  longer  possible 
because  the  aodal  responses  cannot  be  equal  In  aaplltude  but  opposite  in 
sign.  Even  though  node  point  rotating  backward  at  constant  angular 
speed  Is  thus  ruled  out,  It  aay  still  be  possible  to  have  a  node  line 
rotating  at  an  average  angular  speed  of  -w.  The  conditions  required  for 
the  forward  rotating  node  point  are  very  similar  to  the  undamped  case, 
as  will  be  shown. 

Conditions  Required  fqx  g  HadS  gains  •  To  obtain  a 

rotating  node  point  In  the  case  of  the  damped  disk  excited  by  a  rotating 
load  requires  that 


cosa,  -cosaj  Jo\ 

sina,  -sinaj  {Riiui))  \o/ 

Setting  the  determinant  equal  to  zero  yields 

-cosa,  sino2  + Sina,  cosa, -0  -  sin  (a,  -  a,) 

The  only  obvious  solution  to  this  equation  Is 


(127) 


(128) 


a,  -aj 


(129) 


because  the  phase  angles  usually  range  from  0  to  ir  radians,  being  zero 
only  for  «  -  0  and  k  only  for  w  -  •  .  Substituting  this  result  back 
Into  the  first  row  of  matrix  equation  127  gives 


(130) 
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Thus,  for  the  damped  disk  excited  by  a  rotating  load,  the  conditions  for 
a  rotating  node  point  are  that  the  response  amplitudes  be  equal  and  the 
phase  angles  be  equal  as  well. 

There  may  be  other  Instances  In  which  a  node  point  rotating  at 
approximately  the  excitation  frequency  of  the  disk  could  occur,  ks 
noted  In  the  stationary  load  examples,  approximate  satisfaction  of  the 
phase  angle  requirement  may  yield  something  like  a  rotating  node  line. 

In  this  light,  consider  the  case  of  light  damping  or  large  mlstunlng. 

In  such  case,  equation  129  may  not  be  the  only  possible  phase 
requirement  because  It  Is  possible  that 


a,-a3*/r 


(131) 


This  may  well  be  the  more  Important  solution,  since  equation  123  might 
not  be  satisfied  at  a  frequency  where  the  mlstuned  modes  are  dominant. 

As  with  the  case  of  the  stationary  load.  It  Is  desirable  to 
determine  what  the  rotating  node  point  requirements  are  in  terms  of 
modal  damping  values  and  excitation  frequencies.  These  relations  are 
determined  In  the  next  two  subsections. 

Requirement  211  Modal  Damning.  To  solve  for  the  requirements 
on  modal  fractions  of  critical  damping,  equation  121  Is  used  to  generate 
the  following  Identities: 


( 1 +(2Ci/*cu/co,]  -  - -  <^32) 

^  w  V  I  ■/  j  sina,  cosa, 


Substituting  this  equation  into  the  second  row  of  equation  127  gives 


Qor,/,(r,)/,(ro)_^_2_ 

- - r-sin  a 

I  (2^ ,  n(Jo/(xj  1 ) 


‘7or,/2(r,)/2(roj__ 
/  sin 


0 


(133) 
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Then,  eppIlcaClon  of  eqiietion  129  yields 


it, (2^1 /ICO /to  i)  A:2(2f  j/ito/toj) 


(134) 


Solving  this  equation  for  the  fractions  of  critical  damping  yields 


Ca  /a(r,)/2(ro)Ac,co2  ^^35^ 

?i  /i(r,)/,(ro)Ac2CO, 

which  Is  apparently  Independent  of  the  excitation  frequency  (u . 

Requirement  on  Frequency .  Substitution  of  equation  132  Into 
the  first  row  of  equation  127  rather  than  the  second  yields  the 
excitation  frequency  requirement  for  a  rotating  node  point.  It  Is 


CO 


/'■<« 


(136) 


which  Is  Identical  to  equation  113  and  Is  Independent  of  the  modal 
damping  values. 

As  In  the  case  of  the  undamped  disk,  there  Is  no  guarantee  that 
this  equation  will  yield  a  real  frequency,  since  the  response  aaq>lltudes 
need  not  be  equal  In  the  frequency  regime  where  It  Is  possible  for  the 
phase  angles  to  be  equal.  An  example  problem  will  help  to  Illustrate 
these  requirements. 

Rxampi A  Problem  #7 .  Consider  the  sane  parameters  as  In 
example  #6.  Applying  equation  136  obviously  yields  the  same  frequency 
as  equation  115  did  In  example  #6  because  the  equations  are  Identical. 
The  excitation  frequency  Is 

co„,  =  2,093.415  radlans/sec 
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Th«  addition  of  daaping,  however,  introduces  an  additional  requirement 
in  the  form  of  equation  135  .  Applying  this  equation  yields 

^-0.9877 

S  I 

To  plot  the  response,  one  fraction  of  critical  damping  must  be  selected. 
Therefore,  let  -  0.004,  then  f 2  “  0.00395  .  The  response  amplitudes 
and  the  difference  between  them  are  plotted  in  figure  24  in  the 
frequency  band  where  the  rotating  node  point  requirements  are  met. 


^rgsp<ms»  (RAOIANS/SEC) 


Figure  24.  Example  Problem  #7,  Response  Amplitudes 

-  0.004,  fj  -  0.00395;  Forward  Rotating  Node  (?) 
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Ic  nust  be  stressed  Chat  a  rotating  node  point  in  this  case  Is 
quite  unlikely  due  Co  Che  inaccuracy  of  equation  35  at  this  frequency. 

A  more  interesting  case  may  be  one  in  trhich  equation  131  is  employed. 
Such  a  condition  is  considered  as  a  final  example. 

Example  Problem  #£.  In  example  #4  it  was  shown  that 
rotating  node  lines  could  occur  even  if  the  phase  requirement  was  not 
met.  The  result  was  that  the  speed  of  the  mode  was  not  constant  but  had 
an  average  speed  eqiial  to  <i»/n  .  Similar  results  are  expected  in  the 
rotating  load  case.  Consider  the  same  parameters  as  in  example  #1 
except  for  the  torsional  spring  constant.  Let  K  -  50.0  instead  of  24.0. 
This  will  introduce  some  additional  mistiming  so  that  equation  131  can 
be  invoked.  The  resulting  mode  shapes  are  plotted  in  figure  25. 


Figure  25.  Example  Problem  #8,  Mistuned  Radial  Mode  Shapes 
a:  -  50 . 0  :  a/b  -  0 . 1  ;  Mode :  m  -  0 ,  n  -  2 
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The  nodal  fractions  of  critical  danplng  will  still  have  to  be  very 
snail  to  satisfy  equation  131.  The  approach  should  be  Che  sane  as  Chat 
used  In  example  #4.  That  Is,  both  nodal  danplng  values  should  be 
specified,  and  the  response  amplitudes  should  be  sec  equal  for  some 
radial  coordinate  as  follows:  let  -  0.0010,  -  0.0011  and 

Tg  -  O.OS  meters.  Solving  equations  130  and  121  yields  the  following 
excitation  frequency  and  phase  angles: 

cu  =  830.691  rad/sec  a ,  =  3 .03357  a,  »  0 . 1  1 897  J^adlans 

These  conditions  are  suitable  for  equation  131. 

The  node  line  shape  is  plotted  at  Its  extreme  curvature  and 
magnified  100  tines  In  figure  26 .  The  node  line  has  more  curvature  in 
this  case,  as  expected,  because  of  the  sClffer  torsional  springs. 


Figure  26.  Example  Problem  #8,  Node  Line  Shape  Magnified  100  Times 
Mode:  m-0,  n-2  ;  JC-50.0 

As  in  example  #4 ,  the  disk  response  may  be  plotted  at  various 
times  to  determine  whether  or  not  it  has  the  shape  of  the  tuned  mode 
shape  and  hoe  the  displacement  amplitude  varies  as  the  node  line  travels 
around  the  disk.  This  is  done  in  figures  27-29. 

L 


I 
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Figure  27.  Example  Problem  #8,  Disk  Response  Vs.  Theta 

C  1.017  milliseconds;  -  0.001,  "■  0.0011 


Figure  28.  Example  Problem  #8,  Disk  Response  Vs.  Theta 

c  “  1.490  milliseconds;  -  0.001,  {‘2  “  0.0011 
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Figure  29.  Example  Problem  #8,  Disk  Response  Vs.  Theta 

C  -  2.435  milliseconds;  -  0.001,  {'2  ”  0.0011 

These  figures  show  that  the  response  does  have  the  shape  of  the  tuned 
mode  changing  in  amplitude  as  it  rotates,  Just  as  the  disk  of  example 
#4. 

Figure  30  plots  the  node  line  position  vs.  time.  This  figure 
shows  that  the  node  line  speed  variations  are  mild  compared  to  those  of 
example  #4.  Note  also  that  this  is  a  backward  rotating  node  line,  as 
predicted  for  the  undamped  disk.  Thus,  the  case  of  small  damping  is  a 
good  approximation  to  the  undamped  case.  If  additional  damping  were 
introduced,  as  in  a  more  realistic  case,  the  node  line  speed  variations 
could  be  kept  minimal  by  increasing  the  amount  of  mistuning.  This  would 
also  further  distort  the  node  line  shape,  but  it  is  obviously  not 
visibly  distorted  in  figures  27-29. 
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Figure  30.  Example  Problem  #8,  Hode  Line  Foslclon  Vs.  Time 
-  0.001,  -  0.0011;  Hode:  m  -  0,  n  -  2 

The  rotating  load  yields  results  that  are  only  slightly  different 
from  those  of  the  stationary  load.  Differences  in  the  phase 
requirements  for  a  rotating  node  point  result  from  the  fact  that  the 
rotating  load  excites  the  two  modes  out  of  phase.  The  modal  phase 
angles  must  differ  by  0  or  x  with  a  rotating  load.  These  phase  angles 
with  respect  to  the  load  allow  the  two  modes  to  respond  x/2  radians 
apart,  which  is  the  same  relationship  between  the  modal  responses  as 
that  required  for  a  rotating  node  point  in  the  case  of  a  stationary 
load.  This,  however,  virtually  eliminates  the  possibility  of  a  forward 
rotating  node  point  in  favor  of  the  backward  rotating  node  point.  Also, 
the  magnitude  of  the  angular  speed  of  the  rotating  node  lines  is  u 
rather  than  w/n,  as  in  the  case  of  the  stationary  load.  All  other 
aspectj  of  the  rotating  load  response  are  the  same  as  those  of  the 
stationary  load. 
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Th®  rotating  load  Is  also  a  model  for  a  disk  rotating  past  a 
stationary  load.  The  reasons  and  impact  of  using  this  model  fcr  2 
rotating  disk  are  examined  in  the  next  chapter. 
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V.  RotatlnE 


The  equation  of  notion  of  the  rotating  disk  has  not  been  solved  In 
closed  fom,  but  enough  Is  known  about  It  to  estimate  the  forced 
response  with  a  fair  degree  of  certainty.  Proven  estimates  of  the 
natural  frequencies  abound  In  the  literature  (5:1143).  The  angular 
dependences  of  the  mode  shapes  are  generally  assumed  to  be  sinusoidal. 
This  dependence  must  certainly  be  slnxisoldal  In  the  limit  of  small 
rotation  speeds,  since  It  must  converge  to  the  stationary  disk 
eigenfunction  In  this  limit.  Making  use  of  this  assumption  will  allow  a 
development  of  the  mlstuned  rotating  disk  forced  response  that  Is 
directly  parallel  to  that  of  the  mlstuned  stationary  disk.  Because  the 
exact  response  cannot  be  determined,  the  response  will  be  stated  only  In 
a  qxialltatlve  fashion. 

£.aaat.l9n  Motion 

The  eqtiatlon  of  motion  for  the  transverse  vibration  of  a  rotating 
disk  can  be  determined  by  the  methods  used  In  chapter  II.  Only  the 
additional  terms  due  to  the  rotation  need  be  derived  here.  For  the 
rotating  disk,  the  kinetic  energy  of  the  vibrations  and  work  done  by  a 
point  load  are  identical  to  those  of  the  stationary  disk.  The  kinetic 
energy  due  to  disk  rotation  Is  not  considered  because  it  Is  assumed  to 
be  Independent  of  the  disk  transverse  velocity.  The  difference  comes  In 
the  potential  energy  of  the  disk.  The  additional  potential  energy  due 
to  the  rotation  is 


(137) 


El 

where  the  radial  and  angular  stresses  respectively  due  to  a  rotation 
speed  0  are  determined  by  applying  the  center-clamped  boundary 
PI  conditions  to  the  general  stress  equations  given  by  Saada  (10:336). 

They  are 
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(13») 

and 

where 

b^(3^v)-^a"(1^3v) 

and 

r  g'fa'[feV3^y)^a^(I^3v)] 

Substituting  the  energy  term  of  equation  137  Into  Hamilton's  Principal 
and  Integrating  by  parts  yields  three  additional  terms  for  the  equation 
of  motion  and  two  additional  terms  for  the  boundary  conditions.  The 
complete  equation  of  motion  for  the  rotating  disk  is 
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The  boundary  condition  equations  6,  7,  and  10  are  unchanged.  The 
boundary  condition  equations  8  and  9  become 
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The  solution  of  equation  140  is  unknown;  however,  early  (4:276) 
and  recent  (5:1143)  researchers  alike  have  assumed  it  is  of  the  form 


u;(r,0,()  =  a(O/(r)sinrt0  (1^3) 

which  is  the  same  as  the  form  of  the  solution  of  the  stationary  disk 
equation  of  motion.  This  is  a  fair  assumption.  Because  the  two 
solutions  must  converge  in  the  limit  of  small  rotation  speeds,  they  are 
bound  to  have  distinct  similarities.  This  assumption  is  more  than  is 
necessary  to  make  the  mistuning  model  work  for  the  rotating  disk, 
because  the  two  mis tuned  modes  will  decouple  as  long  as  the  integral  of 
the  two  $  dependences  from  zero  to  2ir  vanishes.  Although  the  exact 
natural  frequencies  and  the  radial  dependences  of  the  mode  shapes  are 
unknown,  the  natural  frequencies  can  be  approximated.  This  thesis  will 
not  treat  these  approximations,  as  the  qualitative  behavior  of  the  disk 
is  the  main  focus  here  (see  5:1143  and  11:2010). 

Rotating  Disk  Mistuning 

The  assumption  of  equation  143  allows  the  rotating  disk  to  be 
mistuned  in  the  same  fashion  as  the  stationary  disk.  In  the  case  of  the 


rotating  disk,  tha  foundation  of  the  torsional  springs  nust  rotate  with 
the  disk.  Application  of  the  torsional  spring  alstuning  model  again 
results  in  two  distinct  self*adjoint  eigenvalue  problems  and  two 
distinct  orthogonal  modes  (sine  and  cosine  modes).  The  eigenvalues 
cannot  be  solved  for  in  closed  form,  but  the  forced  response  of  the  disk 
can  still  be  determined  qualitatively. 

Response  £2  d  Stationary  Point  Lsd^ 

Once  sxibjected  to  the  assumption  of  equation  143,  the  forced 
response  of  the  rotating  disk  to  a  stationary  point  load  is  exactly 
parallel  to  that  of  the  stationary  disk  to  a  rotating  load.  Thus,  all 
of  the  equations  and  concepts  of  the  chapter  IV  rotating  load  section 
are  applicable  to  the  rotating  disk.  The  only  major  difference  between 
the  two  responses  is  that  the  rotating  disk  modal  parameters  are  all 
functions  of  the  disk  rotation  speed  and  they  cannot  be  solved  for  in 
closed  form. 

The  differential  equations  of  motion  for  the  rotating  disk  are 
precisely  equations  114  and  115  except  that  the  modal  masses, 
stiffnesses,  and  natural  frequencies  are  now  unknown  functions  of  the 
disk  rotation  speed.  These  equations  are  obtained  by  assuming  that  the 
response  can  be  approximated  by  the  response  of  only  the  two  mistimed 
modes,  as  in  equation  35  (recall  that  this  assumption  is  valid  for  the 
self-adjoint  system  providing  there  are  no  other  resonant  frequencies 
close  to  the  mistuned  modes  being  studied) .  Equation  35  is  substituted 
into  the  equation  of  motion  (eq.  140),  and  the  resulting  equation  is 
multiplied  through  by  each  mistuned  mode.  The  two  equations  resulting 
from  this  multiplication  are  then  integrated  over  the  domain  of  the  disk 
to  yield  two  inde;.endenc  equations.  The  modal  damping  is  then  added  to 
these  equations  to  yield  the  form  in  equations  114  and  115.  Once  the 
parameters  in  these  equations  are  determined  the  remaining  equations  of 
chapter  IV  are  applicable  with  the  load  rotation  speed  0  replaced  by  the 
disk  rotation  speed  (■). 
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The  rotating  disk  response  to  a  stationary  load  has  all  of  the 
features  of  the  stationary  disk  response  to  a  rotating  load.  For  light 
daaplng  and  significant  nistuning,  node  lines  rotating  backward  through 
the  disk  would  be  observed  when  the  nodal  response  aaplitudes  are  equal. 
The  node  lines  of  such  a  response  would  not  be  able  to  nove  through  the 
disk  at  the  constant  angular  speed  Q  but  would  instead  aove  faster  and 
then  slower  with  an  average  angular  speed  0.  Node  lines  rotating 
forward  (stationary  in  the  non- rotating  reference  frame)  in  the  disk 
could  only  be  observed  if  the  nistuned  nodes  were  the  dominant  modes  at 
the  frequency  where  -  Oj,  which  is  highly  unlikely  since  this 
frequency  is  generally  much  higher  Chan  the  natural  frequencies  of  the 
mlstuned  modes. 

The  response  of  Che  rotating  mlstuned  disk  as  given  here  in 
qualitative  form  provides  Che  designer  with  a  reasonable  understanding 
of  Che  resonant  behavior  of  such  a  disk.  Only  the  stationary  point  load 
has  been  addressed,  and  bladed  disks  have  not  been  addressed;  however, 
this  thesis  should  be  useful  as  a  qualitative  tool. 
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VI.  Results  ADji  Discussion 


Mlstunlng  Model 

The  mlstunlng  model  used  In  this  thesis  was  contrived  with 
mathematical  simplicity  as  a  primary  criterion;  however,  It  is  fairly 
flexible  and  allows  for  the  de- coupling  of  the  modal  equations  of 
motion.  In  addition,  the  model  can  be  tuned  to  decrease  the 
Imperfection  and  Its  effects  until  they  are  arbitrarily  small,  which  Is 
a  property  one  would  expect  a  real  disk  to  have. 

In  most  cases,  the  cause  of  mlstunlng  In  a  given  disk  Is  unknown; 
therefore,  no  mlstunlng  model  Is  assured  of  being  correct.  It  Is  also 
unclear  as  to  how  different  mlstunlng  models  affect  the  modal  responses; 
however,  it  Is  possible  to  Imagine  two  models  that  would  have  profoundly 
different  effects  on  the  mathematical  form  of  the  response.  For 
example,  consider  a  model  using  point  masses  Instead  of  the  springs  used 
In  the  present  model.  Such  a  model  would  not  be  self-adjoint,  and, 
therefore,  the  mode  shapes  would  become  series  of  the  eigenfunctions  in 
the  point  mass  system.  The  disk  response  In  the  neighborhood  of  a 
resonance;  however,  may  still  be  similar  to  the  response  of  the  current 
model,  since  a  single  mode  will  dominate  the  others  near  Its  resonant 
frequency.  For  this  reason,  many  reasonable  mlstunlng  models  may  yield 
a  similar  response  In  the  neighborhood  of  a  resonance.  If  this  Is  true, 
then  the  main  difference  among  various  models  may  be  the  complexity  of 
the  mathematical  behavior.  Although  the  model  used  In  this  thesis  does 
not  have  a  form  one  might  expect  in  a  real  disk,  it  may  represent  the 
real  disk  as  well  as  the  next  model,  and  it  has  the  added  benefit  of 
simplicity. 

The  flexibility  of  the  model  is  due  to  the  use  of  variable-stiff¬ 
ness  springs.  Springs  at  the  outer  edge  allow  the  amount  of  mistuning 
to  be  varied  over  a  wide  range,  although  the  mode  shapes  will  differ 
greatly  with  a  large  spring  stiffness,  A  well  manufactured  disk. 
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however,  would  noC  have  enough  mistuning  to  require  such  a  large  spring 
stiffness  in  the  model. 

The  model  is  unique  in  that  it  ensures  orthogonality  between  the 
mistuned  modes.  Because  the  springs  are  placed  at  the  nodes  a  given 
mode,  they  are  not  felt  by  the  ”un-sprung*  mode.  On  the  other  hand, 
since  the  boundary  condition  is  specified  as  free  at  the  nodes  of  the 
"sprung"  mode,  the  "sprung"  mode  cannot  feel  the  free  boundary 
condition.  Note  that  self-adjointedness  of  the  eigenvalue  problem  is 
not  required  to  de-couple  the  differential  equations  of  motion  but  is 
needed  to  ensure  that  an  individual  eigenfunction  can  respond 
independently  of  the  other  eigenfunctions.  Thus,  the  placement  of  the 
imperfections  (springs)  is  what  helps  this  model  work. 

Stationary  Disk  Response 

The  response  of  the  stationary  mistuned  disk  will  usually  have 
rotating  node  lines  in  the  neighborhood  of  a  resonance  and  can  have  the 
approximate  shape  of  the  tuned  mode  shape  under  more  precise  conditions. 
It  was  found  that,  as  long  as  both  mistuned  modes  are  not  in  phase  with 
each  other,  the  node  lines  of  the  disk  must  move.  They  must  move  from 
the  node  lines  of  one  mode  to  the  node  lines  of  the  other,  which  are  at 
the  lines  of  peak  displacement  of  the  first,  and  they  must  do  so  twice 
per  cycle  of  excitation.  This  is  true  whether  the  applied  load  is  a 
stationary  oscillating  point  load  or  a  rotating  point  load.  If  the  two 
modal  response  amplitudes  are  equal,  the  disk  will  take  on  the  shape  of 
a  single  tuned  mode  shape  rotating  around  the  disk;  however,  the 
amplitude  of  this  disk  shape  will  not,  in  general,  remain  constant  in 
time.  Only  if  the  two  modes  respond  exactly  90  degrees  out  of  phase 
will  the  rotating  disk  shape  amplitude  be  constant  in  time.  In 
addition,  the  node  lines  will  move  at  the  constant  speed  d9/dt=  ±w/n  or 
±w  (corresponding  to  the  stationary  and  rotating  load  cases 
respectively)  only  if  both  of  the  rotating  node  point  requirements 
(equal  response  amplitudes  and  90  degree  phase  shift)  are  met. 
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In  general,  both  of  the  rotating  node  point  requirements  cannot  be 
satisfied  simultaneously  for  all  radial  coordinates  of  the  disk,  but,  if 
the  mistuning  is  small  and  one  coordinate  satisfies  the  requirements , 
the  remaining  coordinates  will  almost  satisfy  the  requirements.  In  such 
case,  the  node  line  will  distort  very  slightly  as  it  travels  around  the 
disk,  first  bending  one  way  and  then  the  other.  At  times  when  one  mode 
is  at  zero  response,  the  node  line  will  be  straight.  As  the  node  line 
moves  from  this  time,  it  will  bend  according  to  the  difference  in  the 
response  amplitudes  at  a  given  radial  coordinate.  This  bending  will 
increase  to  a  maximum  when  the  node  line  reaches  a  point  between  the  two 
node  lines  of  the  mistuned  modes,  then  it  will  decrease  to  zero  as  it 
approaches  the  node  line  of  the  other  mode.  On  passing  the  mode  node 
line,  it  will  begin  to  bend  again. 

The  rotating  node  point  requirements  as  stated  in  the  previous 
paragraph  are  obviously  the  same  for  both  stationary  and  rotating  point 
loads,  but  the  two  load  types- have  differences  in  their  phase  shift 
capabilities.  The  phase  shift  between  the  two  modes  is  always  90 
degrees.  However,  if  the  phase  angles  between  the  load  and  the 
responses  are  considered,  the  requirements  appear  to  be  different.  This 
is  because  the  rotating  load  automatically  excites  the  modes  out  of 
phase;  therefore,  for  the  rotating  load,  an  equal  phase  shift  of  each 
mode  with  respect  to  the  load  causes  a  90  degree  phase  shift  between  the 
two  modes.  The  difference  in  the  phase  shift  capabilities  of  the  two 
load  types  results  in  node  lines  that  rotate  in  different  directions. 
Whereas  both  forward  and  backward  rotating  node  lines  are  theoretically 
possible  with  both  rotating  and  stationary  loads,  only  the  backward 
rotating  node  line  is  likely  to  occur  with  a  rotating  load 
(corresponding  to  phase  differences  between  the  modal  responses  of  plus 
and  minus  90  degrees).  The  forward  rotating  node  line  is  unlikely  with 
a  rotating  load  because  the  excitation  frequency  needed  to  obtain  the 
proper  phase  shift  is  usually  too  far  away  from  the  resonance. 
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Rotating  &ia]s  Response 

The  disk  rotating  past  a  stationary  point  load  was  found  to 
respond  much  like  the  stationary  disk  with  a  rotating  load,  If  Its  mode 
shape  Is  assumed  to  vary  sinusoidally  with  $ .  The  solution  of  the 
rotating  disk  equation  of  motion  Is  currently  unknown.  Therefore,  some 
assumption  about  Its  form  must  be  assumed.  The  assumption  that  the  8 
dependence  Is  sinusoidal  has  been  made  by  most  previous  Investigators 
and  seems  reasonable  because  It  converges  nicely  to  the  stationary  disk 
solution  In  the  limit  of  zero  disk  rotation  speed.  Using  this 
assumption  allows  a  qvialltatlve  determination  of  the  rotating  disk 
forced  response.  The  same  mlstunlng  model  Is  used  to  mlstune  the 
rotating  disk,  except  that  the  torsional  spring  foundation  must  move 
with  the  disk.  Once  the  modal  parameters  m^,  k^,  and  are  determined, 
the  response  of  the  disk  rotating  past  a  stationary  load  can  be  found 
using  the  stationary  disk  rotating  load  equations.  The  modal  parameters 
as  well  as  the  radial  variation  of  the  mode  shape  cannot  be  determined 
In  closed  form,  but  approximate  methods  of  evaluating  the  natural 
frequencies  have  been  published  (5:1143,  11:2010),  Thus,  at  least  a 
qualitative  understanding  of  the  rotating  disk  response  Is  Immediately 
available. 


I 


83 


VII.  Conclusions 


By  the  use  of  a  new  model  for  the  forced  response  of  a  mistuned 
disk,  It  was  shown  that  a  mistuned  disk  may  respond  with  rotating  node 
lines  much  like  the  the  response  postulated  by  Tobias  and  Arnold  for  the 
perfect  disk.  A  model  was  devised  to  mlstune  the  disk  without  coupling 
the  eigenfunctions  Into  Infinite  series  of  Bessel  functions.  Then  the 
forced  response  of  the  stationary  disk  to  both  rotating  point  loads  and 
stationary  point  loads  oscillating  In  time  was  determined  and  found  to 
have  rotating  node  lines  under  certain  conditions.  The  disk  rotating 
past  a  stationary  load  was  also  considered  and  found  to  have  a  response 
similar  to  the  stationary  disk  excited  by  a  rotating  load. 

The  new  mlstunlng  model  was  contrived  to  prevent  the  mistuned 
eigenfunctions  from  coupling  between  orders  of  the  Bessel  functions.  It 
uses  torsional  springs  placed  along  the  outer  boundary  of  a 
center -clamped  annular  disk.  One  spring  Is  placed  at  diametrical  node 
of  a  given  mode.  The  boundary  conditions  split  the  mode  Into  two  modes 
that  are  orthogonal  to  each  other  and  to  all  other  modes.  For  the 
stationary  disk,  the  modes  are  each  fully  determined  by  four 
coefficients  of  Bessel  functions  and  Modified  Bessel  functions  of  the 
first  and  second  kinds  all  of  the  same  order. 

The  disk  response  to  a  point  load  can  take  the  form  of  rotating 
node  lines  and,  under  certain  conditions,  may  even  take  the  shape  of  the 
tuned  mode  shape  rotating  through  the  disk.  The  node  lines  must  rotate 
through  the  disk  as  long  as  the  two  modes  respond  out  of  phase  with  each 
other.  If  the  response  amplitudes  are  equal,  the  response  will  look 
like  a  tuned  mode  shape  rotating  through  the  disk.  If,  In  addition,  the 
modes  are  responding  90  degrees  out  of  phase,  the  amplitude  of  the  disk 
shape  will  be  constant  in  time  and  the  node  line  angular  speed  will  be 
precisely  ±w/n  or  ±w  (corresponding  to  the  stationary  and  rotating  load 
cases  respectively).  To  meet  this  last  condition,  the  modal  damping 


muse  meet  stringent  requirements.  These  results  are  also  applicable  to 
the  disk  rotating  past  a  stationary  point  load  If  the  rotating  disk 
eigenfunctions  are  assumed  to  have  a  sinusoidal  9  dependence.  This 
assumption  was  made  because  the  rotating  disk  equation  of  motion  has  not 
been  solved. 
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